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ABSTRACT 

We compute the one loop graviton contribution to the self-energy of a very 
light fermion on a locally de Sitter background. This result can be used 
to study the effect that a small mass has on the propagation of fermions 
through the sea of infrared gravitons generated by inflation. We employ di- 
mensional regularization and obtain a fully renormalized result by absorbing 
all divergences with BPHZ counterterms. An interesting technical aspect of 
this computation is the need for two noninvariant counterterms owing to the 
breaking of de Sitter invariance by our gauge condition. 
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1 Introduction 

In this paper we compute and renormalize the one loop quantum grvaitational 
corrections to the self-energy of very light fermions on a locally de Sitter 
background. The physical motivation for this exercise is to facilitate a later 
study of how inflationary gravitons affect fermions and, in particular, the 
contrast between the case of exactly massless fermions and those with a 
small mass. Nonzero mass introduces two competing effects: it changes how 
fermions propagate and it also alters how they interact with gravity. The 
first of these changes tends to suppress the effects of inflationary gravitons 
because it makes the fermion wave function oscillate so that interactions at 
different times tend to cancel. However, the new interaction enhances the 
effect of inflationary gravitons because it does not fall off with time. 

The current work can be seen as complementing two previous studies of 
massless fermions on de Sitter. In both cases the technique was to compute 
the one loop fermion self-energy — i[jEj](rc; x') and then use it to solve the 
quantum-corrected Dirac equation for fermion mode functions, 



The first model results from Yukawa coupling the fermion to a massless, min- 
imally coupled (MMC) scalar on a nondynamical de Sitter background [3]. 
The second model consists of the fermion with dynamical gravity on de Sitter 
background [1]. Powers of the inflationary scale factor a = e are crucial 
for understanding the results in both cases. The self-energy from the (pipifj 
interaction of the first model grows like a ln(a) relative to the classical term. 
The result fermion mode functions behave as if they had a growing mass. 
The interactions of the second model all possess derivatives — for example, 
dhipif) — which limit the induced self-energy to grow no faster than In (a) 
relative to the classical term. The resulting fermion mode functions behave 
as if they had a growing field strength, which could be understood as the 
random walk that fermions take under buffeting from the sea of inflationary 
gravitons [1, 2]. Although the effect from gravitons is smaller than that from 
massless scalars, it is universal, independent of assumptions about the ex- 
istence or couplings of unnaturally light scalars. It is even conceivable that 
the graviton effect might, in a more complicated model, lead to baryogenesis 
during inflation. 

What we expect for massive fermions in dynamical gravity is that the 
absence of derivatives in the mhipip interaction will cause the self-energy to 
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grow like a ln(a) relative to the classical kinetic term, and like m(a) relative 
to the classical mass term. When the classical mass is large (relative to the 
Hubble parameter) we expect at most a small enhancement of the fermion 
field strength. When the classical mass is small, classical dynamics are mostly 
controlled by the kinetic term and we expect the quantum correction to 
have a much larger proportional impact. One might intuitively expect the 
crossover to come for fermion masses near the Hubble parameter. However, 
we shall specialize to the case of very light fermions, both because this is 
where the largest effects should occur, and because expanding in the fermion 
mass makes an enormous simplification in the computation. 

This work also deserves a place in the growing list of studies of quantum 
infrared effects during inflation. Among these are: 

• The effects of self-interacting, MMC scalars on nondynamical de Sitter 
background [45, 46, 47, 48]; 

• The effects of a charged, MMC scalar on nondynamical de Sitter back- 
ground [81, 83]; 

• The effects of a nonlinear sigma model on nondynamical de Sitter back- 
ground [51, 52]; 

• The effects of a MMC scalar on gravitons on de Sitter background 
[49, 50]; 

• The effects of gravitons on a MMC scalar on de Sitter background [43]; 
and 

• The effects of gravitons on interacting conformal matter on de Sitter 
background [53, 54]. 

It should also be noted that the series of leading infrared logarithms can 
be summed for scalar potential models using the stochastic technique of 
Starobinksy and Yokoyama. The same resummation can be achieved for 
Yukawa theory, and for scalar QED, but it has so far not been accomplished 
for either the nonlinear sigma model, or for quantum gravity. Each fully 
renormalized quantum gravitational result is an important piece of "data" 
in the search for such a resummation. 

Dirac + Einstein is not perturbatively renormalizable [5], however, ultra- 
violet divergences can always be absorbed in the BPHZ sense [6, 7, 8, 9]. A 
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widespread misconception exists that no valid quantum predictions can be 
extracted from such an exercise. This is false: while nonrenormalizability 
does preclude being able to compute everything, that not the same thing as 
being able to compute nothing. The problem with a nonrenormalizable the- 
ory is that no physical principle fixes the finite parts of the escalating series 
of BPHZ counterterms needed to absorb ultraviolet divergences, order-by- 
order in perturbation theory. Hence any prediction of the theory that can be 
changed by adjusting the finite parts of these counterterms is essentially ar- 
bitrary. However, loops of massless particles make nonlocal contributions to 
the effective action that can never be affected by local counterterms. These 
nonlocal contributions typically dominate the infrared. Further, they cannot 
be affected by whatever modification of ultraviolet physics ultimately results 
in a completely consistent formalism. As long as the eventual fix introduces 
no new massless particles, and does not disturb the low energy couplings 
of the existing ones, the far infrared predictions of a BPHZ-renormalized 
quantum theory will agree with those of its fully consistent descendant. 

It is worthwhile to review the vast body of distinguished work that has 
exploited this fact. The oldest example is the solution of the infrared prob- 
lem in quantum electrodynamics by Bloch and Nordsieck [10], long before 
that theory's renormalizability was suspected. Weinberg [11] was able to 
achieve a similar resolution for quantum gravity with zero cosmological con- 
stant. The same principle was at work in the Fermi theory computation 
of the long range force due to loops of massless neutrinos by Feinberg and 
Sucher [12, 13]. Matter which is not supersymmetric generates nonrenormal- 
izable corrections to the graviton propagator at one loop, but this did not 
prevent the computation of photon, massless neutrino and massless, confor- 
mally coupled scalar loop corrections to the long range gravitational force 
[14, 15, 16, 17]. More recently, Donoghue [18, 19] has touched off a minor 
industry [20, 21, 22, 51, 52] by applying the principles of low energy effective 
field theory to compute graviton corrections to the long range gravitational 
force. Our analysis exploits the power of low energy effective field theory in 
the same way, differing from the previous examples only in the detail that 
our background geometry is locally de Sitter rather than flat. 1 

That summarizes why the exercise we have undertaken is both valid and 
interesting. In the next section we work out the fermionic part of the Feyn- 
man rules for Dirac + Einstein. Section 3 is devoted to the issues associated 

1 For another recent example in a nontrivial cosmology see [25]. 
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with the graviton propagator. A major complication concerns the impossibil- 
ity of employing an average, de Sitter invariant gauge condition [26, 27]. We 
give a short review of the complex literature on this issue. Then we introduce 
a noninvariant gauge fixing term, isolate the subgroup of de Sitter transfor- 
mations that it respects, and present the gauge-fixed graviton propagator. 
The section closes with a discussion of the BPHZ counterterms necessary for 
our computation. In section 4 we evaluate the contributions from diagrams 
involving a single 4-point interaction. In section 5 we evaluate the more dif- 
ficult contributions which involve two 3-point interactions. Renormalization 
is accomplished in section 6, and our conclusions are given in section 7. 

2 Fermions in Quantum Gravity 

The coupling of gravity to particles with half integer spin is usually accom- 
plished by shifting the fundamental gravitational field variable from the met- 
ric 9nv{ x ) to the vierbein e^ m (x). 2 Greek letters stand for coordinate indices 
and Latin letters denote Lorentz indices, and both kinds of indices take val- 
ues in the set {0, 1, 2, ... , (D — 1)}. One recovers the metric by contracting 
two vierbeins into the Lorentz metric i] bc , 

g^u(x) = e IJ , b (x)e uc (x)r] bc . (2) 

The coordinate index is raised and lowered with the metric (e% = g^e^), 
while the Lorentz index is raised and lowered with the Lorentz metric (e M b = 
V bCe fj,c)- We employ the usual metric-compatible and vierbein compatible 
connections, 

e/3fe ;M = =>> A licd = e v c (e vd>li -r f ' lu/ e pd ) . (4) 

Fermions also require gamma matrices, 7^. The anti-commutation rela- 
tions, 

{ 7 6 , 7 C }^( 7 V + 7 C 7 6 ) = -2^ C /, (5) 

imply that only fully anti-symmetric products of gamma matrices are actually 
independent. The Dirac Lorentz representation matrices are such an anti- 

2 For another approach see [28]. 
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symmetric product, 



j* c = l( 7 y _ 7 y) = l 7 iy . (6) 



4V ' ' ' / 2 



They can be combined with the spin connection (4) to form the Dirac co- 
variant derivative operator, 

V, = d, + l -A, cd J cd . (7) 

Other identities we shall often employ involve anti-symmetric products, 

7 6 7 y = 7 ty 7 d l _ ri bc ^ d + ^y - ri cd ^ b ? (8) 
^jcd = l 7 [y 7 d] + l^y _ l^y _ ^ 

We shall also encounter cases in which one gamma matrix is contracted into 
another through some other combination of gamma matrices, 

7%, = -DI, (10) 
7 y 7b = (£>-2) 7 c , (11) 
7Y7S = 4r7 cd 7 - (L»-4) 7 V , (12) 

7 y 7 y 7b = 2 7 y 7 c + (D-4) 7 y 7 e . (13) 

The Lagrangian of massive fermions is, 

C m ^ c = i)e\^ b W^^g -m^^g . (14) 

Because our locally de Sitter background is conformally flat it is useful to 
rescale the vierbein by an arbitrary function of spacetime a(x), 

ep b = aep b e pb = a" 1 e pb . (15) 

Of course this implies a rescaled metric g^ v , 

9,u = a 2 g, u g^ = a~ 2 r u - (16) 

The old connections can be expressed as follows in terms of the ones formed 
from the rescaled fields, 

r% = ^'K^+n^-r^vj + fV (17) 

A^d = -a~ 1 (^ c e lxd -e" d e IM ?)a jl/ + A^a . (18) 
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Although the massive Dirac is not conformal invariant, we still re-scale 
fermion fields to make it simpler, 

* = a^i) and ^ = a^i) . (19) 
After performing the above rescalings, the Dirac Lagrangian is 

£ D irac = ^z\l h iV^s[-g- ypij , (20) 

where = <9 M + f A^ cd J cd . 

One could follow early computations about flat space background [29, 30] 
in defining the graviton field as a first order perturbation of the (conformally 
rescaled) vierbein. However, so much of gravity involves the vierbein only 
through the metric that it is simpler to instead take the graviton field to be 
a first order perturbation of the conformally rescaled metric, 

g pu = rj^y + Kh^ with k 2 = IQttG . (21) 

We then impose symmetric gauge (e^ = e^) to fix the local Lorentz gauge 
freedom, and solve for the vierbein in terms of the graviton, 

e[#W = (vW 1 )/ = Vf}b + \ Kh ? b ~ ^ h p h ib + ■■■ (22) 

It can be shown that the local Lorentz ghosts decouple in this gauge and one 
can treat the model, at least perturbatively, as if the fundamental variable 
were the metric and the only symmetry were diffeomorphism invariance [31]. 
At this stage there is no more point in distinguishing between Latin letters for 
local Lorentz indices and Greek letters for vector indices. Other conventions 
are that graviton indices are raised and lowered with the Lorentz metric 
(h^ = rf p h pv , h pv = rf p r\ V(J h pa ) and that the trace of the graviton field is 
h = r] pu h pv . We also employ the usual Dirac "slash" notation, 

fij = V^j • (23) 

It is straightforward to expand all familiar operators in powers of the 
graviton field, 

e\ = 6% - \nh\ + ?« V V + ■ ■ ■ , (24) 

gl» = rf» _ ^ + K ^ p h PV , (25) 

\[-g = 1 + \kh + ln 2 h 2 - \n 2 h pcJ h pa + . . . (26) 
2 8 4' 
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Applying these identities to the conformally rescaled Dirac Lagrangian gives, 

K- 



£Dirac — ^ 
+ 



tqj — am 
1 



2 



hiyl-h^^idv-h^Y J pa - om/i 



-h 2 --h pa h 



pa 



-f///i i/i 2 - ^ Q/3 /i Q J**l + C(/t 3 ). (27) 
8 4 ^ J 

We observe that the rescaled massive fermion propagator in a locally de 
Sitter background could be derived from the solution of Candelas and Raine 
[32, 33]. The relation between them is only up to some powers of scale factors, 

iS[m}(x;x') c .R. = (aa'y^iSlmKx;^), (28) 
therefore the conformally re-scaled fermion propagator is, 



iS[m](x] x) 



H D ~ 2 ^ { D x. ,. D -i, .^1 fa A 



(4tt 



r(f-i+^)r(f-^ 



m\"p/D • m "> 



X 



m D ,m D y\/I~l ( 



I r(f-i)r(f) 
r(^-i-^)r(#+^ 



„ /D id u in u ii 

FA l-i— , — h«—; — ; 1-- 



if 2 if' 2 

m .D ,m D y\/i+7 



r(f-i)r(f; 



if' 2 if' 2 



■■ (29) 



where the coordinate interval is Ax 2 (x;x') = — of' || 2 — ( j 77 — 77' j — iS) 2 and 
%Tp here is just a~(^~\g)a(^~\ 

It is useful to recast the solution (29) using the transformation formula 
of hypergeometric functions [34] and expands out its series expansion, 



iS[m](x; x 
4 



r(f-i: 

A 2. 



iq) + am 



Ax D-2 

(fiW)f^r(f-i)r(2-f)(if) 



(47r )f r(i+*f)r(i-^) 



^+(^-l)^ffa 7 ° + 



am 



x 



n=0 1 



r(n+f-l+»g)r(n+f-l-»g) 

r(n+f-i)r(n+i) 

nn+l+i^nn+l-i^U i% 



' H 



+ 7 L 



(I)' 



W\n+2- 



Notice that (30) recovers the fermion propagator in flat spacetime when 
taking H = 0, 



lim i 

h — >0 



iSj 



[x; x) = i 



iSj 



cf 



[x; x') + % 



iSj 



fin 



X', x ) , 



iSj 



iSj 



r(^-i) 1 



fin 



>; a:') 



r(f-l) ma 



4tt 2 



(31) 
(32) 

(33) 



Here "cf " stands for "conformal" and "fm" stands for "flat spacetime mass" . 

Because we only endow fermions very small mass compared with the 
Hubble parameter, for the computation purpose we simplify the infinite series 
expansion by only keeping terms at order m, 



iSj 



(x; x') 



iSj 



cf 



(x; x') + % 



iSj 



fm 



(x; x') 



oo 



//W)^r(f-i)r ( 3-f) %0 + ( |_ i)/fa 



H> (4tt) 
T(n- 



D 
2 



1) /■yw 



n=0 1 



(!) 



(2-f) 

T(n + 1) 



r(n + l) V4^ r(n + 3-f)M 



+ C(m 2 ). (34) 



Even though the two infinite series expansions tend to cancel out with each 
other in D = 4, the combinations are still finite owing to the divergent 
factor l D . In addition, they can not be expressed in closed form such as 

Hankel function. The facts which we mentioned above complicate the whole 
computation. 

We now represent the various interaction terms in (27) as vertex operators 
acting on the fields. At order k the interactions involve fields, tyj and 
h a p, which we number "1", "2" and "3", respectively. Each of the three 
interactions can be written as some combination V"^ of tensors, spinors and 
a derivative operator acting on these fields. For example, the first interaction 
is, 



Hence the 3-point vertex operators are, 



(35) 



K 



V 3ij 



0: 



VI 



3/x 



K 

——amrj 



<xf5j 



IJ • 



(36) 
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-\ K 2 rf^idp 2 
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— ^K 2 r) a ^r) p(T am 


5 




10 


^K?r} ap rfP am 



Table 1: Vertex operators Uj(j pa contracted into \E r i\J/ jh a ph i 



The order k 2 interactions define 4-point vertex operators Ujfj Pa similarly, for 
example, 

\n 2 h 2 m^ = ^VV^, x ^jh^h^ = Uff a x H^l-h^h^ . (37) 

o o 

The ten 4-point vertex operators are given in Table 1. Note that we do not 
bother to symmetrize upon the identical graviton fields. 



3 Graviton Propagator and Counterterms 

The gravitational Lagrangian of low energy effective field theory is, 

^Einstein = ~[^q{ R ~ {D-2)h)y/=jj . (38) 

The symbols G and A stand for Newton's constant and the cosmological 
constant, respectively. The unfamiliar factor of D — 2 multiplying A makes 
the pure gravity field equations imply R^ v = Ag^ u in any dimension. The 
symbol R stands for the Ricci scalar where our metric is spacelike and our 
curvature convention is, 

r = g^R^ = gr (rv, - r v + r vr V - r? va r° m ) . (39) 

Even though gravity is not conformally invariant, it is still useful to express 
it in terms of the rescaled metric (16) and connection (17), 

^Einstein = { a°- 2 R- 2 (D - l) a D " 3 ^ (a^ - f 

-(D-4)(D-l)a D - 4 r u a»a u -(D-2)Aa D ^^g~ . (40) 
9 



The factors of a which complicate this expression are the ultimate reason 
there is interesting physics in this model! 

We need to fix a in order to work out the graviton propagator from the 
Einstein Lagrangian (40). The unique, maximally symmetric solution for 
positive A is known as de Sitter space. In order to regard this as a paradigm 
for inflation we work on a portion of the full de Sitter manifold known as the 
open conformal coordinate patch. The invariant element for this is, 

ds 2 = a 2 (—di] 2 + dx-dx) where a(rj) = ——— , (41) 



and the D-dimensional Hubble constant is H = \JA/(D — 1). Note that the 
conformal time f] runs from — oo to zero. For this choice of scale factor we 
can extract a surface term from the invariant Lagrangian and write it in the 
form [27], 

Einstein -Surface = (f-DF^-^fg^V,/.^ + ^ ' spgT^^ 
x | ^hapjjhu^^ \h a p^ph a p^ v + \h a p php Va \h a p^php cr ^ . (42) 

Gauge fixing is accomplished as usual by adding a gauge fixing term 
directly to the Lagrangian. One can add such a gauge fixing term and then 
use the well-known formalism of Allen and Jacobson [35] to solve for a fully 
de Sitter invariant propagator [36, 26, 37, 38, 39]. However, a curious thing 
happens when one uses the imaginary part of any such propagator to infer 
what ought to be the retarded Green's function of classical general relativity 
on a de Sitter background. The resulting Green's function gives a divergent 
response for a point mass which also fails to obey the linearized invariant 
Einstein equation [26]! We stress that the various propagators really do 
solve the gauge-fixed, linearized equations with a point source. It is the 
physics which is wrong, not the math. Furthermore, quantum corrections 
bring new problems when using de Sitter invariant gauges. The one loop 
scalar self-mass-squared has recently been computed in two different gauges 
for scalar quantum electrodynamics [44]. When the computation was done 
in the de Sitter invariant analogue of Feynman gauge the result was on-shell 
singularities! Off shell one-particle-irreducible functions need not agree in 
different gauges [55] but they should agree on shell [56]. 

The nature of the problem is the apparent inconsistency between de Sitter 
invariance and the manifold's linearization instability. Any propagator gives 
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the response (with a certain boundary condition) to a single point source. If 
the propagator is also de Sitter invariant then this response must be valid 
throughout the full de Sitter manifold. But the linearization instability pre- 
cludes solving the invariant field equations for a single point source on the 
full manifold! This feature of the invariant theory is lost when a de Sitter in- 
variant gauge fixing term is simply added to the action 3 . To understand this 
subtle issue, we digress to explain the difference between an "exact gauge" 
and an "average gauge" . The former is obtained by choosing the gauge pa- 
rameter to make the fields obey some equation at each point in space and 
time. The other more common type of gauge fixing which we have already 
mentioned in the preceding paragraph is accomplished by adding some terms 
to the invariant Lagrangian. The functional integral representation for this 
type of gauge condition can be viewed as a weighted "average" of exact 
gauge. One can usually find the equivalence between the exact gauge condi- 
tions of the canonical formalism and the average gauge condition employed 
in the functional formalism. The surprising thing is that the Faddeev-Popov 
ansatz proven by the Sidney Coleman [57] for the flat space on the manifold 
R A is not generally valid for the gauge theory which possesses a linearization 
instability [58]. 

Because of this topological obstacle to adding covariant, average gauge 
fixing terms on gravity in de Sitter, except for employing the exact gauge 
condition on an open manifold [59] one can also avoid it either by working 
on the full manifold with a noncovariant gauge condition that preserves the 
elliptic character of the constraint equations; or else by employing a covariant, 
but not de Sitter invariant gauge on an open submanifold [27]. We choose 
the later course and employ the following analogue of the de Donder gauge 
fixing term of flat space, 

Cgf = - V-y^F, , F, = 7/" (V* - \ h P^ + (D-mah.X)- (43) 

Because our gauge condition breaks de Sitter invariance it will be neces- 
sary to contemplate noninvariant counterterms. It is therefore appropriate 
to digress at this point with a description of the various de Sitter symme- 
tries and their effect upon (43). In our D-dimensional conformal coordinate 
system the \D(D + 1) de Sitter transformations take the following form: 

3 It results in changing the constrained equation from being elliptical to hyperbolic, 
which also alters the causal structure of the de Sitter geometry. 
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1. Spatial translations — comprising (D — l) transformations. 

rf = V, (44) 
x H = x* + e\ (45) 

2. Rotations — comprising \{D — l)(D — 2) transformations. 

rf = 77, (46) 
a* = . (47) 

3. Dilatation — comprising 1 transformation. 

77' = krj, (48) 
x H = kx 1 . (49) 

4. Spatial special conformal transformations — comprising (D—l) trans- 
formations. 

rf = , (50) 

i-2e-x+e 2 x-x 



x H = \ * . (51) 

l-26-x+\\0\\ 2 x-x 



It is easy to check that our gauge condition respects all of these but the 
spatial special conformal transformations. We will see that the other sym- 
metries impose important restrictions upon the BPHZ counterterms which 
are allowed. 

It is now time to solve for the graviton propagator. Because its space and 
time components are treated differently in our coordinate system and gauge 
it is useful to have an expression for the purely spatial parts of the Lorentz 
metric and the Kronecker delta, 

V = ^ + <W and C = ^-«- (52) 

The quadratic part of ^Einstein + C-gf can be partially integrated to take the 
form ^h^D J^hpa, where the kinetic operator is, 

D pa — i - fi ct) - -n n pa - 1 S p S a 1 D a 

u p,v — \ 2 °n °v ^V/j-uV 2(£> — 3) M J 



-WW Db + - { ^ ) D C , (53) 
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and the three scalar differential operators are, 

d a = ^(V^ru 



1 /L>-2s 



Do - 



(54) 
(55) 

(56) 



The graviton propagator in this gauge takes the form of a sum of constant 
index factors times scalar propagators, 



pv>Ap a 



(57) 



I=A,B,C 

The three scalar propagators invert the various scalar kinetic operators, 

Di x iAi(x; x) = i5 D (x - x) for I — A,B,C , (58) 

and we will presently give explicit expressions for them. The index factors 
are, 



rpA 

V-v 1 pa 

rpB 

rpC 
pa 



— 2 T lp{p r la)v ~ 



D-3 



VauV 



pv I per 



2 



(59) 
(60) 



(£>-2)(£>-3) 



+ V (^-3)W + V -( 61 ) 



With these definitions and equation (58) for the scalar propagators it is 
straightforward to verify that the graviton propagator (57) indeed inverts 
the gauge-fixed kinetic operator, 



D 



pa 



pv 



X I 



pa^ 



[x; x')=5^5^i5 D (x-x') 



(62) 



The scalar propagators can be expressed in terms of the following function 
of the invariant length £(x; x') between x^ and 

y(x;x') ee 4 sin 2 {^H£(x; x')) = aa'H 2 Ax 2 (x;x) , (63) 
= aa'H 2 (\\x-x'\\ 2 -(\r)-r)'\-i5) 2 ) • (64) 
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The most singular term for each case is the propagator for a massless, con- 
formally coupled scalar [60], 

«M*^ = gJr(f-i)(l)*-\ (65) 

The A-type propagator obeys the same equation as that of a massless, min- 
imally coupled scalar. It has long been known that no de Sitter invariant 
solution exists [61]. If one elects to break de Sitter invariance while preserv- 
ing homogeneity (44-45) and isotropy (46-47) — this is known as the "E(3)" 
vacuum [62] — the minimal solution is [45, 46], 

x') = iA c f (x; x') 




Note that this solution breaks dilatation invariance (48-49) in addition to 
the spatial special conformal invariance (50-51) broken by the gauge condi- 
tion. By convoluting naive de Sitter transformations with the compensating 
diffeomorphisms necessary to restore our gauge condition (43) one can show 
that the breaking of dilatation invariance is physical whereas the apparent 
breaking of spatial special conformal invariance is a gauge artifact [63]. 

The B-type and C-type propagators possess de Sitter invariant (and also 
unique) solutions, 



iAsix; x) = iA c f(x; x) 



iA c {x; x) = iA ci (x; x) + 




(67) 



2 

D 



-( n -J + 3 ) r(n+2) (4) [ (68) 

These expressions might seem daunting but they are actually simple to use 
because the infinite sums vanish in D = 4, and each term in these sums 
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goes like a positive power of y(x] x'). This means the infinite sums can only 
contribute when multiplied by a divergent term, and even then only a small 
number of terms can contribute. Note also that the 5-type and C-type 
propagators agree with the conformal propagator in D = 4. 

In view of the subtle problems associated with the graviton propagator 
in what seemed to be perfectly valid, de Sitter invariant gauges [26, 27], it is 
well to review the extensive checks that have been made on the consistency 
of this noninvariant propagator. This propagator has been shown to obey 
the Ward identity at tree order [65] and one loop [72]. In addition, the only 
fully dimensionally regulated loop results for quantum gravity on de Sitter 
background have been obtained using it [1, 43, 66, 67] as well. 

It remains to deal with the local counterterms we must add, order-by- 
order in perturbation theory, to absorb divergences in the sense of BPHZ 
renormalization. The particular counterterms which renormalize the fermion 
self-energy must obviously involve a single ip and a single ip. At one loop order 
the superficial degree of divergence of quantum gravitational contributions 
to the fermion self-energy is three, so the necessary counterterms can involve 
zero, one, two or three derivatives. These derivatives can either act upon 
the fermi fields or upon the metric, in which case they must be organized 
into curvatures or derivatives of curvatures. We will first exhaust the possible 
invariant counterterms for a general renormalized fermion mass and a general 
background geometry, and then specialize to the case of the lowest order of 
mass term in de Sitter background. We close with a discussion of possible 
noninvariant counterterms. 

All one loop corrections from quantum gravity must carry a factor of k 2 ~ 
mass -2 . There will be additional dimensions associated with derivatives and 
with the various fields, and the balance must be struck using the renormalized 
fermion mass, m. Hence the only invariant counterterm with no derivatives 
has the form, 



With one derivative we can always partially integrate to act upon the ip field, 
so the only invariant counterterm is, 



Two derivatives can either act upon the fermions or else on the metric to 
produce curvatures. We can organize the various possibilities as follows, 




(69) 




(70) 





(71) 
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Three derivatives can be all acted on the fermions, or one on the fermions 
and two in the form of curvatures, or there can be a differentiated curvature, 



K\ m (R^-^g^R)^h m W v iPy/^ , K^R^^V^g ■ (72) 

The counterterms with zero mass (72) have already been discussed in [1]. 
Because we are dealing with the order m contribution, we do not consider 
counterterms (69) and (70) and only focus on (71) for the purpose of the 
current computation. The specialization of the invariant counter-Lagrangian 
we require to de Sitter background at order m is therefore, 

A£ inv = X^rmpiiT^f^j^g + \ 2 K 2 mR^^g , (73) 
— > \ lK 2 Vm(i9a-H?)v + \ 2 k 2 {D - l)DH 2 ma^ . (74) 

Here Ai and A 2 are D-dependent constants which are dimensionless. The 
associated vertex operators are, 

C Vl3 = X 1 K 2 [—d 2 + mH-f°^)=\ 1 K 2 m^a~ 1 ^, (75) 
C 2ij = X 2 D(D -l)K 2 H 2 ma . (76) 

Of course C 1 is the higher derivative counterterm mentioned in section I. It 
will renormalize the most singular terms — coming from the iA c f part of the 
graviton propagator — which are unimportant because they are suppressed 
by powers of the scale factor. The other vertex operator, C 2 , is a sort of 
dimensionful field strength renormalization in de Sitter background. It will 
renormalize the less singular contributions which derive physically from in- 
flationary particle production. 

We do not use the background field technique in background field gauge 
[75, 76, 77, 78]. We must therefore come to terms with the possibility that 
divergences may arise which require noninvariant counterterms. What form 
can these counterterms take? Applying the BPHZ theorem [6, 7, 8, 9] to the 
gauge-fixed theory in de Sitter background implies that the relevant coun- 
terterms must still consist of k 2 times a spinor differential operator with the 
dimension of mass-cubed, involving no more than three derivatives and act- 
ing between \& and \&. As the only dimensionful constant in our problem, 
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powers of H must be used to make up whatever dimensions are not supplied 
by derivatives. 

Because dimensional regularization respects diffeomorphism invariance, it 
is only the gauge fixing term (43) that permits noninvariant counterterms. 4 
Conversely, noninvariant counterterms must respect the residual symmetries 
of the gauge condition. Homogeneity (44-45) implies that the spinor dif- 
ferential operator cannot depend upon the spatial coordinate x l . Similarly, 
isotropy (46-47) requires that any spatial derivative operators di must either 
be contracted into 7* or another spatial derivative. Owing to the identity, 

( 7 ^) 2 = -V 2 , (77) 

we can think of all spatial derivatives as contracted into 7*. Although the 
temporal derivative is not required to be multiplied by 7 we lose nothing by 
doing so provided additional dependence upon 7 is allowed. 

The final residual symmetry is dilatation invariance (48-49). It has the 
crucial consequence that derivative operators can only appear in the form 
a _1 9 M . In addition the entire counterterm must have an overall factor of a, 
and there can be no other dependence upon rj. So the most general order m 
counterterm consistent with our gauge condition takes the form, 

A£ non = K 2 H 2 ma$>S({Ha)- 1 ^d G , (tfa)" 1 ^)* , (78) 

where the spinor function S(b, c) is at most a second order polynomial func- 
tion of its arguments , and it may involve 7 in an arbitrary way. 

Three more principles constrain the order m noninvariant counterterms. 
The first of these principles is that the fermion self-energy at order m involves 
only even powers of gamma matrices. This follows because the three-point 
vertices, the four-point vertices and the fermion propagator all consist of an 
even number of 7's at order m 1 and an odd number of 7's at order m°. The 
diagram which consists of one 4-point vertex possesses an even number of 

4 One might think that the they could come as well from the fact that the vacuum 
breaks de Sitter invariance, but symmetries broken by the vacuum do not introduce new 
counterterms [79]. Highly relevant, explicit examples are provided by recent computations 
for a masslcss, minimally coupled scalar with a quartic self-interaction in the same locally 
de Sitter background used here. The vacuum in this theory also breaks de Sitter invariance 
but noninvariant counterterms fail to arise even at two loop order in either the expectation 
value of the stress tensor [45, 46] or the self-mass-squared [47]. It is also relevant that 
the one loop vacuum polarization from (massless, minimally coupled) scalar quantum 
electrodynamics is free of noninvariant counterterms in the same background [80]. 
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gamma matrices at order m. The contribution from any diagram with two 3- 
point vertices consists of three factors involving gamma matrices: one factor 
from the fermion propagator and one factor from each of the two vertices. 
At order m such a product consists of one even and two odd factors, so 
it contains an even number of gamma matrices. This principle fixes the 
dependence upon 7 and allows us to express the spinor differential operator 
in terms of just ten constants fa, 



K 2 H 2 maS{ {Ha)- l -i°d , (if a) "Yd, 



= K 2 mal faia'Vdo) 2 + fa (a" YO))(a~ Y$) + faipT^di) 



+tf 7 U ftfa'VU) + &(a-V^) + H 2 fa . (79) 



In this expansion, but for the rest of this section only, we define noncommut- 
ing factors within square brackets to be symmetrically ordered, for example, 

(a-VdoXa-V^)] = i(a- 1 7°9 )(a- 1 Y^) + ^aT^d^fd^ . (80) 

The second principle is that our gauge condition (43) becomes Poincare 
invariant in the flat space limit of H — > 0, where the conformal time is 
i] = —e~ Ht /H with t held fixed. In that limit only the three quadratic terms 
of (79) survive, 



lim K 2 H 2 maS((Ha)- 1 ^d , {Ha)- 1 1 i d l ) 



n 2 ma{ Ma-Vdo) 2 + fa (crV3))(a~V0i) + faia'^dif \ (81) 



Because the entire theory is Poincare invariant in that limit, these three 
terms must sum to a term proportional to (7 M <9 M ) 2 , which implies, 



fa = 2& = fa ■ 
But in that case the three quadratic terms sum to give (75), 



(82) 



K 2 ma{ (a-^do) 2 + 2 (crYdb)(a 7*$) + (a - V^) 2 = K 2 m$a~ l $ . (83) 
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Because it is the same as one of the invariant counterterms, it need not be 
included in S. Besides, the final term in (79) recovers the other invariant 
counterterm (76). So the two remaining noninvariant counterterms we need 
to consider in (79) are, 

ACnon = ^|/3 4 (a- 1 7 9 ) + ft(a-y3i)}* . (84) 

However, these two terms are not independent of the last term in (75). There- 
fore we could chose any four independent counterterm operators we need for 
this computation, 

777- 

a^n 2 — d 2 , a^K 2 H 2 ma , (85) 
a 

a 2 K 2 mHdo , a^mH^ $ . (86) 



4 Contributions from the 4-Point Vertices 

In this section we evaluate the contributions from 4-point vertex operators 
of Table 1. The generic diagram topology is depicted in Fig. 1. The analytic 
form is, 



■T 



4pt 



10 

i=i 



Iij 



(x; x) 5 (x—x') . 



(87) 




x 



Fig. 1: Contribution from 4-point vertices. 

And the generic contraction for each of the vertex operators in Table 1 is 
given in Table 2. 

From an examination of the generic contractions in Table 2 it is apparent 
that we must work out how the three index factors [apT^] which make up the 
graviton propagator contract into rj al3 and i] ap . For the A- type and 5-type 
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T 
1 


l [al3 L ^pa\\X ) X ) IU j \X~X ) 


l 


— i[ a A p J(x; x) ip o [x — x ) 


2 


^/c z[ ap A aj a](x; x) <pb \x — x) 


3 


1 2 'Try A 1/ \ /? OtT £"/}/ A 

«[ a Q ,A pcr J(x; x) 7 p a d> (x— x J 


4 


— g/c i[~gA a(T J(x; xj 7 p a o (x — x J 


5 


-\k 2 d' p i[ a a A pa ](x; x') yj ap S D (x-x') 


6 


I* 2 ^i[ a /3 A Q(T ](x; x') YJ P ° S D (x-x') 


7 


{k 2 a M i[ a /3 A aa ](x; x) 7" J*" 5 D (x-x') 


8 


\k 2 &H[ afi A pu \{x; x') YJ aa 5 D (x-x') 


9 


-|/t 2 ami[ a a A p p ](x; x) 5 D (x-x') 


10 


jK 2 ami[ al3 A a/3 ](x; x) 5 D (x—x') 



Table 2: Generic 4-point contractions 



index factors the various contractions give, 



7/ 



afS 



rpA 

a/3 1 pa 



(^-3) v 



ap 



a/3 1 pa 



= 



V 



ap 



rpA 

aft 1 pa 



a/3 1 pa 



( D -^)^*> (88) 



= -p-l)« + 5^,(89) 



For the C-type index factor they are, 



77 



rj 



ap 



rnC 

ocP 1 pa 



rpC 

aft 1 pa 



= (-*-) 5050 + 4 



(D-2)(D-3) /pCT ' 
> 2 
" 2 (^) W+ (jD _2)( jD -3)^ 



•D-2- 



(90) 



At order m we actually only require double contractions. For the A-type 
index factor these are, 



^afi^pa 



r) ap rf a 



rpA 

aft 1 pa 



afiT 



(91) 



The double contractions of the 5-type and C-type index factors are, 



^aff^pa 



rpB 

a/3 1 pa 







^ap^fia 



rpB 

a/3 1 pa 



2(D-1), 



(92) 
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I 


J 


i[ Q/3 T/J iAj(x- x') iUf p ° S D (x-x') 


9 


A 


|/t 2 pz|) « m iA^ix; x) S D (x—x') 


9 


C 




10 


A 


\k 2 [D(D 1) 2g~Jj] am2A A (i;i)J fl (i x') 


10 


B 


i/c 2 (D - l)amiA s (a;;x) 5 D (z-x') 


10 


C 


f k2 (d-2) 5 (d-3) am ^A c (x; x) 5 D (x-x') 



Table 3: 4-point contribution from each part of the graviton propagator at 
order m. The vertices 1-8 could only give the contribution at order m°. 



rpC 

ap 1 pa 



(D-2)(D-3) ' 



rpC 

aft 1 pa 



2 P 2 ~5^+8) 

-^( jD _2)( jD -3) • ^ 



Table 3 was generated from Table 2 by expanding the graviton propagator 
in terms of index factors, 



afiApa 



(x; x') 



rpA 

ap 1 pa 



iA A (x; x') + aP T* iA B (x; x') + 



Tc 
.. . , 



pa 



iA c (x;x') . 
(94) 

We then perform the relevant contractions using the previous identities. 

From Table 3 it is apparent that we require the coincidence limits on each 
of the scalar propagators. For the A-type propagator these are, 



lim iA A (x; x) 

rpl \rr> 



H D - 2 T(D - 1) 



f r( 



-7T 



cot(|l>) +21n(a)| . (95) 



The analogous coincidence limits for the 5-type propagator are actually finite 
in D = 4 dimensions, 



lim iA B (x;x') 



H D-2 pp _ x) 



• D N 



X 



D-2 



(4tt)t m 2 

The same is true for the coincidence limits of the C-type propagator, 



(96) 



lim iAc(x; x') 



h u - 2 r d - 1 



X 



( 47 r)T T(f) (D-2)(D-3) 



(97) 
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I 


J 


i am 5 D (x—x') 


9 


A 


(D-l) A 


9 


B 





9 


C 


1 

(D-2) 2 (D-3) 2 


10 


A 


-\[D{D-1)-^]A 


10 


B 


1 (D-l) 

2 (D-2) 


10 


C 


1 (D 2 -5D+8) 

2 (D-2) 2 (D-3) 2 



Table 4: The 4-point contributions at order m. All contributions are multi- 
plied by K2H V T i D ul ] . Here A = % cot(^) -ln(a). 
V J Uk\% r(f) 2 v 2 ) \ i 



We apply the various coincidence limits to each contraction in Table 3 and 
present the order m, 4-point contributions in Table 4. The total summation 
for this local contributions is quite simple, 



S 4pt 



iK 2 H 2 ma H 



D-4 



2 ^+%f (D-4) r(f) 

12 In a- 1 



S D (x-x') 



ik 2 H 2 
+ ; - - ma 



16tt 2 



[x — X 



(98) 



5 Contributions from the 3-Point Vertices 

In this section we evaluate the contributions from two 3-point vertex opera- 
tors. The generic diagram topology is depicted in Fig. 2. The analytic form 
is, 



i3pt 



x; x 



a/3. 
ik 1 



X) I 



,x 



i=i 



;*')E^V) 
j=i 



(99) 
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iVf(x) i[S](x; x') iVj P °(x>) i[ a pA pa ](x; x') 


l 


1 




l 


2 


\K 2 & p {$i[s\(x- x') r*i a a \*}& x 1 )} 


l 


3 


{iK 2 @i[S](x;x')yj a »d' ti i 




l 


4 


^iK 2 am$i[S](x; x') i[ a a A p p ](x; x') 


2 


1 


l^d'^d? i[S}{x; x')Yi 


a pA<> p ){x-x')} 


2 


2 




2 


3 


-\m 2 7 Q <9^[S](x; x') YJ^^iUpAp^x; x') 


2 


4 


— ^iK 2 am yd 13 i[S](x; x') i[ a pA p ](x; x') 


3 


1 


-\iK 2 d' u {-i a JP» i[S](x; x') 7"^ 1 




](x;x')} 


3 


2 


\in 2 d' p {y i[S] (x; x') 7 a a M i 


aflApfj 


(x;x 7 )} 


3 


3 


-i/t 2 7" Jto i[S}(x; x') yj^d^fy i 


apApa (x, X ) 


3 


4 


— \K 2 am y i\S\{x\ x') d^i 




(x; x') 


4 


1 


— ^iK 2 d' p {am i[S](x; x') 7 M i 


[ a a A p p ](x;x')} 


4 


2 


\iK 2 d' p {ami[S](x]x')Y % 


a a A pa }(x;x')} 


4 


3 


— \k 2 am i[S](x] x') yj Gp dp i 






4 


4 


-\K 2 a 2 m 2 i[S}{x;x')i[ a a A p p 





Table 5: Generic Contributions from the 3-Point Vertices. 




Fig. 2: Contribution from two 3-point vertices. 

Because there are four 3-point vertex operators in (36), there are sixteen 
vertex products in (99). We label each contribution by the numbers on its 
vertex pair, for example, 



I-J 



= iV, 



a/3, 



X X I 



S 



[x-x') x iV p /{x') x i 



apApa- 



(100) 
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Table 5 gives the generic reductions, before decomposing the graviton propa- 
gator 5 . Most of these reductions are straightforward but one subtlety deserve 
mention, that is, derivatives on external lines must be partially integrated 
back on the entire diagram. This happens whenever the second vertex is 
J — l or J = 2, for example, 



2-2 



IK 



X I 



a/3 A pa 



s 



[x; x') 7 P i 



(x; x') }. 



(x;x'), (101) 
(102) 



Another simplification we might use for later contractions is that the Dirac 
slash of the conformal part of fermion propagator gives a delta function, 



S 



cf 



x; x 



iS D (x - x') 



(103) 



5.1 Conformal Contributions 

The key to achieving a tractable reduction of the diagrams of Fig. 2 is that 
the first term of each of the scalar propagators iAj(x;x') is the conformal 
propagator iA ci (x; x). The sum of the three index factors also gives a simple 
tensor, so it is very efficient to write the graviton propagator in the form, 



pvA pa 



2 

[x;x') = [27] Kp r] a)u - — -7)^7)^ iA ci (x; x') 



I=A,B,C 



i5Aj(x;x'), (104) 



where i5Aj(x; x') = iAj(x; x')—iA ci (x; x'). In this subsection we evaluate the 
contribution to (99) using the 3-point vertex operators (36) and the fermion 
propagator (34) but only the conformal part of the graviton propagator, 



pv A pu 



(x;x) — > \^lr]p {p r] a)v - jj—^r] IIU r] pa \iA ci (x;x') 



rpcf 

aft 1 pa 



iA ci (x;x) . 
(105) 

We carry out the reduction in three stages. In the first stage the conformal 
part (105) of the graviton propagator is substituted into the generic results 



5 We would not consider the 4-4 contraction because it is an order m 2 contribution. 
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from Table 5 and the contractions are performed. We also make use of gamma 
matrix identities such as (9) and, 



Yjto + ryPjPi* = l -(y»n^ + 7 V /1 )-«7V^ , l*J a ^ = -^P-1)7 M - (106) 

We do not at this stage act any derivatives on the fermion propagator. The 
results of these reductions are summarized in Table 6. Because the conformal 
tensor factor [ a /?T^] contains three distinct terms, and because the factors of 
^ajPn j n Table 5 can contribute different terms with a distinct structure, we 
have sometimes broken up the result for a given vertex pair into parts. These 
parts are distinguished in Table 6 and subsequently by subscripts taken from 
the lower case Latin letters. 

In the second stage we substitute the conformal part of the graviton 
propagator, 

^A c^ (x;x') = J^i^, (107) 

and decompose the fermion propagator (34) into the conformal part, the 
flat spacetime mass term, n = part and n > 1 part of the infinite series 
expansion 6 . In the final stage we act the derivatives. We start from the most 
singular contribution in Table 6, which substitutes the conformal parts of the 
fermion propagator into the contraction 1-4, 2-4, 3-4, 4-1, 4-2 and 4-3 7 . The 
contraction 1-4 and 2-4 vanish owing to the equation (103) and owing to the 
zero contribution from D powers of the coordinate separation in Dimensional 
regularization. We also must remember that [S](a;;a: / ) will be used inside an 
integral in the quantum-corrected Dirac equation (1). For that purpose the 
most singular term at x /fl = x M is quadratically divergent in D = 4 dimensions. 
Hence we first conveniently 8 employ the following identities to express the 
rest of them as a less singular form, 

1 (9 2 1 At — <9 1 

' (108) 



Ax 2D ~ 2 2(D - 2) 2 Ax 2D ~ 4 ' Ax 2D ~ 2 2(D-2)Ax 2D 



-4 



6 We will explain why we separate the n = part from the rest of the infinite series 
expansion in a later paragraph. 

7 The contraction 4-4 is an order m 2 contribution. 

8 Some individual term is easier written as a derivative with respect to acting upon 
a less singular coordinate separation than taking the derivative directly. 
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Table 6: Contractions from the iA cf part of the Graviton Propagator. 
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4 (D-2) 


4 


3 
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Table 7: iA cf (x;x') x i[5'] c f(x; x'). All contributions are multiplied by 
3§ 7 r(f)r(f-l)ma(aa0 1 ^. 



The individual result is quoted in Table 7 and collect all terms of this 
class, 

K 2 ma(aa'f-^ D_ 

1 2 ' 1 2 J 

' ^-?(«-^o-^7°^-r-Ll.(109) 




(D-2) 4 v y u 4 ' r \ Ax 2D ~* 

The expression (109) is still logarithmically divergent in D = 4 after pulling 
out various derivatives. To further renormalize this divergence we extract 
derivatives with respect to the coordinate x^ again, which can of course be 
taken outside the integral in (1) to give a less singular integrand, 

^ 2 t4^- (no) 



Ax 2D ^ 2(D-3)(D-4) Ax 2D 

Expression (110) is integrable in four dimensions and we could take D = 4 
except for the explicit factor of 1/(D — 4). Of course that is how ultraviolet 
divergences manifest in dimensional regularization. We can segregate the 
divergence on a local term by employing a simple representation for a delta 
function, 



d 2 , 1 x _ d 2 f 1 /i D - 4 \ iAnTfj, -* 5 1 

VAx20-eJ /) 1 1 A^20-6 \,./> 2 I ' 



(x—x 1 ) 



D-4VAx 2D - 6 ^ J D-4\Ax 2D - 6 Ax d ~ 2 \ r(f-l) D-4 
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+ 0(D-A)\ + 



r(f-i) D-A 



111) 



After substituting (110) and (111) into (109) one can get, 
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(D — 1) /i D " 4 J -2 m Q2 _ mH f u 



~\n{\i 2 Ax 2 ) 
7i 

mH 



'■ 2 (ua 



Svrf 2 (-D — 2) 
2 G-^o-^7°?p(a;-x'). (112) 



(D-3) (£> - 4) \ (£> - 2) 

To reach the same expressions as the counterterms we mentioned in the 
suction 3 we make use of the following identities, 

^-d 2 5 D (x - x') = {-d 2 + 2Hd \5 D (x - x') , 

-d 5 D (x - x') = \-d - h\5 d (x - x') , 
a 1 a > 

ln{a')d 5 D {x - x') = {ln{a)d + Ha}5 D {x - x') , 
l -^ld 5 D (x - x 1 ) = {^9 + if (l - ln(a)) }5 D (x - x') , 

l ^ld 2 5 D (x - x') = \^-d 2 + 2/f(ln(a) - lU + H 2 a\5 D (x - x') . 
a' 1 a \ / j 

(113) 

After applying (113) to (112) and expanding out (aa') 1_ ~ we get the total of 
this most singular class which is consistent with our counterterm convention, 



^-icfcf 

1 



x; x') = 



D-4 



iK 2 r(f-i)(D-i)/x- 

16ttT (D-3)(D-A) 1 (D-2) a ~ {D-2) 



—2 m „ 9 4mff _ 
a- — _, <9 



mH 1 °^+-mH 2 a \5 D (x-x')^ 



IK 



8tt 2 



ln(a) 



3 m 
~2a 



d 2 +3mHd 



3 — 3 

H — mH^°& mH 2 a 

4 4 



K 



\l™d 2 + -mH{--l)d + -mH 1 ^\d 



3 3 ] 

— -mH do + -mH 2 a >5%x— x') 

ln^Arr 2 ) 



64tt 2 \ 2 a 



Ax 2 



• (114) 



A less singular contribution comes from the flat spacetime mass term of 
the fermion propagator. Note that the contraction (3-3) involves two deriva- 
tives acting upon the conformal graviton propagator, which would produce 



28 



(I-J)sub 


Qi 1 






Pi tla 


,Qf Ha' 








7 ^ Ax 2D-4 


a Ax 2D - 4 


7 ^Ax»-* 


Aa .2D-4 


(1-1) 


D 


-2D 


-2D 


-D 


D 


-D 


(D-2) 2 


(D-2) 2 


(D-2) 2 


2(D-2) 


2(D-2) 


(D-2) 


(1-2) 


-1 


2 


2 


1 


-1 


1 


(D-2) 2 


(D-2) 2 


(D-2) 2 


2(D-2) 


2(D-2) 


(D-2) 


(1-3) 


— rD— •n 


CD— 1 1 






— CD— 11 




2(D-2) 2 


2(D-2) 2 


2(D-2) 2 


4(D-2) 


4(D-2) 


2(D-2) 


(2-1) 


^ 


2 


2 


"I 
1 


1 

— ± 




(D-2) 2 


(D-2) 2 


(D-2) 2 


2(D-2) 


2(D-2) 


(D-2) 


(2-2) a 


— 1 


1 


— 1 


1 

8 


1 

8 


1 
4 


4(D-2) 


2(D-2) 


2(D-2) 


(2-2) b 


—D 


_D 





D 
"8 





D 

T 


4(D-2) 


2(D-2) 


(2-2) c 


1 


-1 


-i 


-1 


l 


-i 


2(D-2) 2 


(D-2) 2 


(D-2) 2 


4(D-2) 


4(D-2) 


2(D-2) 


(2-3) a 


D 


-D 


U 


— D 
16 


U 


-D 
8 


8(D-2) 


8(D-2) 


(2-3) b 


1 


-1 


-1 


-1 


1 


-1 


4(D-2) 2 


4(D-2) 2 


4(D-2) 2 


8(D-2) 


8(D-2) 


4(D-2) 




1 


-1 


1 


-1 

16 


-1 

16 


-1 

8 


8(D-2) 


8(D-2) 


8(D-2) 


(3-l)+(3-2) b 


(2D-3) 
4(D-2) 2 


(2D-3) 
4(D-2) 


(2D-3) 
4(D-2) 


-(2D-3) 
8(D-2) 


(2D-3) 
8(D-2) 


(2D-3) 
8 


(3-2) a 


-1 


-1 

8 


1 

8 


1 


1 


-1 
16 


8(D-2) 


16(D-2) 


16(D-2) 


(3-2) c 


-D 


— D 

8 





D 
16 





-D(D-2) 


8(D-2) 


16 


(3-3) a 





D 

32 





— D 
32 





D(D-2) 
32 


(3-2) b 





1 

32 


-1 

32 


-1 
32 


-1 

32 


(D-2) 
32 


(3-3) c 





-(2D-3) 
16(£»-2) 


-(2D-3) 
16(D-2) 


(2D-3) 
16(D-2) 


-(2D-3) 
16(D-2) 


-(2D-3) 
16 



Table 8: iA c f x ifiSJfm. Note that all contributions are multiplied by the factor 
T^r(f)r(f -l)ma{aa'f-^. 
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a delta function, 



D N 



(D-2) 



A 2. 

4tx 2 



Hha'AT] 2 H 2 aa 



Ax D 



D_ 

ll-K 2 



2Ax D ~ 2 r 



D ' 



8 D (x-x') 



• (H5) 



This delta function would give zero when it is multiplied by D powers of the 
coordinate separation, which occurs in this case. Next, in order to facilitate 
our computation, we make use of the following identity except for (108), 

H 2 aa'Arf Ha- Ha' 



Ax 2D - 2 " 2(D-2) 9 °Ax 2D - 4 ' 



(116) 



and breaks up Ha(a') ^7° into Ha(a')d and Ha(a') r y° f). The intermediate 
results are summarized in Table 8. When all terms in Table 8 are summed, 
we employ (110) for making the expression integrable in D = 4 and also 
segregate the ultraviolet divergences into the local terms using (111). The 
total result we get from this contribution is, 



—1 



-icffm 



IK 



x: x 



\nH 2r(fV 



D-4 



+ (b A -b 2 )Ha }5 D (x-x')^ 



16ttt (D-3)(D-4) 
iH^mH 



(b 2 +h)d +(b 2a +b 3a h°$ 



64tt 2 



ln(a) 



-7V Ha 



+ \Ha 



{x-x') 



K^nH { 
+ 64vr 4 { 



la 3 

Or, a - 

8 a' ° 32 ° 



.A« 7 o^_lZ + l Ha 
16 a' 1 W 32 ' * 4 



d 2 



ln(/i 2 Ax 2 



+ 



3 1 ^ . Ha# 



8 Ha' 4 

b 2 ,b 2a ,b 3 ,b 3a and 64 are D dimension-dependent coefficients, 

(3D - 4) 



Ax 2 
1 1 



Ax 2 



(117) 



b 2 = 

b 2 a = 
h 
ha 



8(D-2) 2 { ' 
(D-4) | (-3D + 5) | 



8(L>-2) 2 
(4D - 7) 



32(D - 2) 
(QD - 17) 
2(D-2) 2 ' lQ(D-2) 
2D -5) 



3 

32 



32(D - 2) 
(D-4) 



(D-4) + 



16(D -2) ' 
(6D - 7) _ 



32(D-2) 16(D-2)' 



-L> + 6 
32~ 



+ 



8(D-2) 



1 



m _ 4) _ _ . 
v ; 16 



118) 
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[i— J ) sub 


a i 


,,()£» 1 


Ha 


Ha' 


H 2 aa'~f°j k Ax k 




7 9 A x^-4 


Ax 2D-4 


A x 2D-4 




(^-^) 


AD 


AD 
(D-2)' 


n 


n 


n 




-A 
(D-2)2 


-4 
(D-2)' 


n 


n 


n 




-(D-l) 
(D-2)2 


-(D-l) 
(D-2)2 


-(D-l)(D-3) 
(D-2)' 


n 


n 


(2-1 ) 


-4 
(D-2)2 


-4 
(D-2) 2 


U 


U 


U 


(2~2) a 


(D-3) 
2(D-2) 


U 


(2D-9) 
2(D-2) 


U 


U 


(2-2) b 


(D-5) 
2 


U 


D(D-4) 
2(D-2) 


U 


U 


(2-2) c 


2 

(D-2)2 


2 

(D-2)' 











(2-3) a 


4 





— 1 
4 


-(D-2) 
A 





(2-3) b 


1 

2(D-2) 2 


1 

2(D-2)' 


{LI- A) 
2(D-2)' 








(2-3) c 


_1 

"T 





_1 

T 


—3 

~ 


_1 
~4~ 


(3-l)+(3-2) b 


-(2D-3) 


-(2D-3) 


-(D-5)(2D-3) 


-(2D-3) 


(2D-3) 


2(D-2) 


2(D-2) 


4(D-2) 


4(D-2) 


4(D-2) 


(3-2) a 


-1 
4(D-2) 


(D-3) 
4(D-2) 


-(D-3) 
8 


1 

8 


1 

8 


(3-2)c 


D-3 
4 





-D 
8 


D 

8 





(3"3)a 








-(D-2)(D-3) 
16 


(D-2)(D-3) 
16 





(3-3) b 





-1 

8 


(D-2) 
32 


(3D-8) 
32 


D-2 
32 


(3-3) c 





(2D- 3) 
4(D-2) 


(2D-3)(3D-7) 
16(D-2) 


(2D-3) 
16 


-(2D-3) 
16 



Table 9: iA c { x i[S] n= o. Note that all contributions are multiplied by the 
factor ^£r(f )r(f - \){aa'f-^ . 
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At the next step we are going to consider contributions from the infinite 
series expansion of the fermion propagator. Because the series carries at least 
one power of mass we only need to consider diagrams which do not originate 
from the mass term in the Lagrangian. Because the infinite series is vastly 
more complicated than other parts of the fermion propagator it would be 
desirable to carry out the computation in D = 4 dimensions. Whether or 
not it is legitimate for us to do this entirely depends on whether this kind of 
contraction is integrable in four dimensions. The dimensionality of the series 
of the fermion propagator is Ax D--i-in an d the one from the conformal part 
of the graviton propagator is Ax d- 2 ■ Also remember that all the terms in 
Table 5 which derive from two order m° vertices carry two derivatives . This 
means that the total dimensionality in this class is a 2D-3-2n • Therefore we 
shall separate the n = part, which working on an arbitrary D dimension is 
necessary, from the rest of the infinite series expansion, which is integrable 
in four dimensions. Because n = part is not integrable in D = 4 it worths 
mentioning its simplification from (34) by performing ^ expansion for gamma 
functions rather than expanding it out around D — 4, 



r(f- 



.771 ,D 



1) +0{m 2 ), 



m 



l±i^(l) +0{m 2 ). 
ri J 



(119) 
(120) 



Here ^'s stand for digamma functions and they cancel out completely at 
order m when one substitutes the above equations back to the n = part of 
series, 

h d-2 r(f-i) 



s 



y X, X ) n = 



(47r)f (2-f) 



(,_)( aa ) 2 



,,o 



X 



+ 



r(3-f)r ( £. 

1 Ha 
(2-f) Ax D ~ 



f mHaa' 



167T? 



Ax D - 



mH D - 3 (aa') 



(47T) 



-r(?)r(2- 



D 



)Ha 



(121) 



Note that the final two terms which have the same D = 4 limit in (121) 
tend to cancel out with each other. The first derivative of (121) has the same 
pattern, 



d,[S] n -_ 



v 2 



!)• 



Svrf 



7^7 
Ax D - 



(D-2) 



YtAx v Ax^ 
Ax® 



+ 



Ax D - 
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(I - J) sub 




£)2 i?a' 


Ha<y°d @ A J D - (i 




tt2 „2 u,.l 


Had Q Ax 2D-e 


^0 Ax 2D - (i 






Ax 2D ~ 6 


(1-1) 


2D 





2D 





2D 


(D-2) 2 (D-3) 


(D-2) 2 (D-3) 


(D-2) 2 
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_2 





_2 





_2 


(D-2) 2 (D-3) 


(D-2) 2 (D-3) 


(D-2) 2 


(1-3) 















2(D-2) 2 (D-3) 


2(D-2) 2 (D-3) 


(2-1) 


2 





2 





2 


(D-2) 2 (D-3) 


(D-2) 2 (D-3) 


(D-2) 2 


(2-2) a 














2(D-2) 


(2-2) b 


— 1 











— D 


2(D-3) 


2(D-2) 


(2-2) c 












1 


(D-2) 2 (D-3) 


(D-2) 2 (D-3) 


D-2 


(2-3) a 


8(D-3) 














(2-3) b 










\ 





4(D-2) 2 (D-3) 


4(D-2) 2 (D-3) 


(2-3) c 


U 


1 


U 


U 


U 


16(D-3) 


(3-l)+(3-2) b 


-(2D-3) 





-(2D-3) 





-(2D-3) 


4(D-2)(D-3) 


4(D-2)(D-3) 


8(D-2) 


(3-2). 














1 

16 


(3-2) c 


D-2 
8(D-3) 











D 
16 


(3-3) a 


-(D-2) 


(D-2) 











32(D-3) 


32(D-3) 


(3-3) b 





1 





1 





32(D-3) 


32(D-3) 


(3-3) c 


(2D-3) 





(2D- 3) 








16(D-2)(D-3) 


16(D-2)(D-3) 



Table 10: iA c{ x i[S] n=0 . Note that all contributions are multiplied by the 
factor ^r(f)r(f-l)(aa') 2 ^- 
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+ 



HaAx, 



2#W 



Ax D ' 2 (D-4)Ax D ~ 



m(H 2 aa' 



_D 

2 " 



-r(5+i)r( 2 -f)^ 



(122) 



The final two terms of equations (121) and (122) would give a non-zero 
contribution when they are multiplied by the divergent term 9 . 

The results derived from this class are lengthy and we would tabulate 
them separately based on their distinctive characteristics. Some contractions 
would produce at least one (D— 4) factor. One source of (D— 4) is from total 
derivatives acting upon (aa') 2 ~~. This factor can arise when one power 
of (aa') comes from the fermion propagator and the rest of it, (aa') 1- ", 
originates from the conformal part of the graviton propagator, i.e. (1-1), 
(1-2), (2-1), (2-2), (3-1), (3-2). Another source of (D-4) comes from the 
following peculiar gamma function contraction, 



yy 7 °7/3 = (D-4) 7 V + 2(D-2)?7 



0;/ 



(123) 



which occurs in the contractions (2 — 2) b , (2 — 3) a , (3 — 2) c and (3 — 3) a . We 
summarized the terms without (D— 4) 10 in Table 9 and Table 10 whereas the 
terms with the (D— 4) factor are presented in Table 11. 

After segregating the divergences into the local terms, the total result 
from Table 9 and Table 10 is, 



— i 



-icfnO-l 



D-4 



X\X 



m 2 ff 2 2r(|) fj, 
32vrf (D-3) (D-4) 



m r 



ld 1 —^^+(d 2 + d 3 + d 4: )ma 



x 5 u (x-x ) + — - ln(a) 



+ 



32vr 2 

9m(d+a') -| ^ 2 |-ln(/i 2 Ax 



m n ^ 9 
-2— 77+-ma 
H 4 



5\x-x')^ 



k 2 H 2 



64tt 4 



32 



Ax 2 



— 5 m 3 13 . 

On + -ma-\ ma 

16 H 4 32 



d 2 



1 



9m(a+a') „ 9 

■On 



3m(5a+a') _ 



16 °Ax 2 16 ' ~ ur Ax 2 

Here d±, d 2 , d 3 and d 4 are dimension- dependent coefficients, 
(D — 8) (3D — 4) 



(124) 



d 1 



8(D-2) 2 



9 One can consult the various gamma function contractions with (122) in Appendix A. 
10 The contractions (1—3), (2— 3)b, (2— 3) c , (3— 3)b and (3— 3) c produce no (D-4) factor. 
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(I— J) sub 
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(1-1) 








4D 
(D-2)2 








(1-2) 








(D-2Y 
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3 

2(D-2) 


(D-2) 


(2-2) b 


-(D-5) 
2(D-2) 


-(D-5) 
2(D-2) 


-(3D-4) 
2(D-2) 
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(D-3) 
4(D-2) 


(D-3) 
4(D-2) 


1 

4(D-2) 
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8 


1 

4 
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(3-l)+(3-2) b 








-(2D-3) 
4(D-2) 


-(2D-3) 
4(D-2) 


(2D-3) 
4(D-2) 


(3-2) a 








-1 
8 


-(D-4) 
8(D-2) 


1 

8 


(3-2) c 


-(D-3) 
4(D-2) 


-(D-3) 
4(D-2) 


-(D-2) 
4) 


(D-3) 
4 





(3-3) a 








{D-3) 
16 


-(D-3) 
16 


-(D-2) 
16 


(3-3) b 

















(3-3) c 


















Table 11: iA cf x i[S] n=0 with a (D — 4) pre-factor. All contributions are 
multiplied by ^#r(f )r(f - l)(aa') 2 ^^. 
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da = 



4 (D - 2, 
;D-6) r l (2D 



3) 



16 



32 8(D-2)J 



(D-2) 



d A = (D 

16 V 



1). (125) 



Because the contributions from Table 11 carries a factor of (D—4), they would 



survive only if they combine with a 



m-4) 



after making them integrable in 



D = 4 dimensions. The result from this part is finite so we give an expression 
for the D = 4, 

k 2 H 2 



—i 



-icfhO-2 



[x; x 1 ) ■ 



64tt 4 



8 if ' ^ 32 32 I 



Ax 2 



(126) 



The contributions which originate from the final two terms of (121) and 
(122) tend to cancel. They could be taken special care first before they 
are tabulated because most of contributions from this class give finite re- 
sults in D = 4 dimensions except for a few divergent terms from the double 
derivatives on A c f(x;x'). The strategy to deal with the finite part of the 
contributions is to perform the (D—4) expansion and make use of (111) and 
the following key identities, 



d, 



d, 



>*Ax D - 2 
1 



-■fi 



D-4 



d 



"Ax 2D - 6 
1 



2D- 



1 + -(£>— 4)(1 — In Arc 2 ) 

1 



d, 



^ Ax 2 ' 



i4n 2 



r(f-i) 



Ax D-2 

Here we present one example from (3 — 2) c , 



l + -(D-4)(2-21nAa; 2 ) 
S D (x-x'). 



1 



Ax 2 



(127) 
(128) 

(129) 



,D^,D .^m./c 2 / H 2 (aa') 



1- — 

z 2 



r( 2 )r( 2 X){ -h\\32^\(D-4) 



DHa 
2(D-3) 



d 2 + ^H 2 a 2 d { 



H D - 2 D £(3-f) 
2^+%^ l 2 >\(D-4) 



(D^2) d + J H ° do 




There are three distinctive contributions in (130). We extract out the pre- 
factor r(y)r(y— 1) to avoid the repeated and lengthy expressions. The 
first kind comes from the two single derivative terms, 



DH 2 a 2 



\(D-4) 2 



1 + i (D -4) (- In ad + 2 - 2 In fi 2 Ax 2 ) 



+ 



-l + -(D-4) 
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(I— J) sub 




7°^S£ 


H 2 a 2 Ha' 
Ax 2 


Had 2 [<iff 2) ] 


(1-1) 


-2-41n(|) 


-2-41n(f) 
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-I -i MI) 
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8 


(3-2). 


-I MI) 


-I MI) 





1 

8 


(3-2) c 


-IMI) 








1 

2 


Sub— total 





I -MI) 









Table 12: iA ci x i[S , ] n=0 form the two terms which tend to cancel. All 
contributions are multiplied by K 64 m f . 
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(I - J) sub 
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H 2 aa'Ar] 


H 2 a 2 -/°-y k Ax k 


H 2 aa'^°-/ k Ax k 


H 2 a 2 Ha' 
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Ax 4 
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1 
2 


1 

2 








1 
4 


(3-3) b 


1 
8 


1 

8 


1 

8 


1 

8 


1 

16 


(3-3) c 


5 
8 


5 
8 


5 
8 


5 
8 


-5 
16 


Sub— total 








3 
4 


3 
4 






Table 13: iA c{ x ifS'j^o from the two terms which tend to cancel. All 
contributions are multiplied by K ^.™^ x (l + ln(|)). 



H 2 

x(-ln — --l-l + lnu 2 Ax 2 ) 
4fi 2 



DH 2 a 2 



Ax 2 



8 



-In 



da 



Ax 2 
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and we presented the temporal and spatial contributions separately in Ta- 
ble 12. The second kind is the contributions which produce a delta function 
originated from the two double derivative terms, 



DHa \ 
i(D-4) 2 \ 
i4n~5 D (x—x') 



l + i(L>-4)(-bW-2)] + [_l + I(£>-4)(- In ^ 



1 + 1} 



x 



r(f-i) 



-DHa 



In + 1 

2/x 



«7T 2 /i 1 



r(f- 



■i) 



-<5 D (:r — x) 



(132) 



and the final results are displayed in the first column of Table 14. The third 
kind is the residual term from d 2 



l 

Ax 2D- 



-fr r x — (D-4) d 2 

\{D-A) 2 y 1 8 



In n 2 Ax 2 



Ax 2 



(133) 



and we give the results in the final column of Table 12. 

The final contribution in this category is from (3 — 3) which consists of 
some finite terms and local divergent terms. Recall in (115) that two deriva- 
tives acting on the conformal part of the graviton propagator would produce 
a delta function. In general, it would be zero in dimensional regularization 
when it acts on dimension-dependent power of the coordinate separation. 
However, the last term with the divergent coefficient in (121) does not pos- 
sess any dimension-dependent power of Ax 2 . As a result, when i^^o are 
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(I - J) sub 


l ^g-ma5 A (x - x') 


f^m^rnaS^x-x') 


(3-l)+(3-2) b 


5[l + ln(f )] 





(3-2) a 


-[i + Mf)l 





(3-2) c 


-4[l + ln(f)] 





(3-3) a 





r(f ) x f 


(3-3) b 





r(f)x(i-f) 


(3-3) c 





r( D 2 i) x ^ 



Table 14: iA cf x i[S , ] n=0 . The delta function contribution from the two terms 
which tend to cancel. 



multiplied by d.dA c f, no any other terms in the calculation can be used to 
cancel this particular divergent local term. Here we present (3 — 3) a as an 
example, 



H 



D-2 



r(f)r(3-f) 



-{D-2) mK 2 \ H\aa') 2 -T 
~T6 H~2\Z2-k d \{D-A) 

DH 3 a 2 a'Ar] 2 D H 3 a 2 a 



DH 3 a 2 a'Ar] 2 D H 3 a 2 a' 



Ax 2D ~* 



2 Ax 2D ~ e 



2 D + 1 n D \{D-A) 



Ax D 



2 Ax D ~ 2 



DHai2iT2 n 
+ — S D (x-x') 



r(f) 



(134) 



We employ the same trick to deal with the finite part, extract out the pref- 

actor r(#)r(g-i) £ff-<£«> 



64n D 



16 



1 + - (D - 4) (- In aa' - 2 In fi 2 Ax 2 ) 



+ 



-l + l(D-A)(-\n^- 2 -l+\n^Ax 2 ) 



x 



2D-8 



\{D-A) 



DH 3 a 2 a'An 2 D H 3 a 2 a' 
' f 



Ax< 



2 Ax 2 



, (135) 



and tabulate the result for D = 4 in Table 13. The divergent term can be 
read off directly and we present it in the second column of Table 13. Finally 
we enclose this sub-class by summing up all the terms from Tables 12, 13 
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and 14, 



—i 



^cfnO-3 



X] X 



^- 2 r(f)r(f-i)r(3-f)(2D-3) 



IK 



ik 2 H 2 ma -a. .. k 2 H 2 \ 
-d\x—x)-\ 



32tt 2 4 



2(£>-4) 
64tt 4 \ 



ma5 D (x—x') 



8 8 V 



ma 2 



8 



l + ln(^) mtfaa' . (136) 

4 J Aar 



~ 1 



The last computation in this sub-section involves the rest of the infinite 
series expansion which are all integrable and hence we can compute it in 
D = 4 dimensions directly. The fermion propagator for the n > 1 of the 
series in four dimensions we employed is, 

mHaa'^ 



167T 2 

and its derivative 11 is, 



n=l 
41 



Ax 2 



+n ln(f )]+# a[l+(n+l) ln(f )] } (137) 



mHaa' y j ^-f°Ax v Ax^ 



Ax 4 



27^7° 



Ax 2 



1+nln- 

4 



+ 



2//a(^7 J/ 7°Ax l/ +Ax jU ) 



Ax 2 



(2n-l) + (n 2 -n) In 



(2n+l) + (n 2 +n) ln^ 



+# 2 a 2 <5j|(2n+3) + (n 2 + 3n+2)ln|| |> . (138) 



One interesting pattern is that taking the derivative of the coefficient 
of the logarithm term with respect to n gives the coefficient for the term 
without logarithms. Before we table the result from each length expression, 
we present the result from the contraction (1-1), 

2K%mS](x;x')^A c{ (x-x')} h ' 2 '" H 



327T 4 



Ax 6 



Ax 6 



x 



(3n 2 -2n-2) + (n 3 -n 2 -2n) ln(^) 



2# 2 a 2 7 VA^ 



AHa 

Ar 4 



(3n 2 -l) + (n 3 -n) ln( 



Ax 4 



(3n 2 +4n-l) + (n 3 +2n 2 -n-2) ln(^-) 



+ 



H 2 a 2 Ha' 



Ax 4 



(3n 2 +6n+2) + (n 3 +n 2 +2n) ln(- 



*4 J 



4' 
(139) 



L One can find various gamma functions contracted with (138) in Appendix A 
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(I-J)sub 


y ~f k Ax k 
Ax 6 


Ari 
Ax G 


Ai 6 


Ai 6 


tfcxAr? 2 
Ax 6 


Ha'Ar, 2 
Ax 6 


(1-1) 


16(3n 2 -2n-2) 


-16(3n 2 -2n-2) 


16(3n 2 -2n-2) 





-16(3n 2 -2n-2) 





(1-2) 


-4(3n 2 -2n-2) 


4(3n 2 -2n-2) 


-4(3rt 2 -2rt-2) 





4(3n 2 -2n-2) 





(1-3) 


6(2n+l) 


-6(2n+l) 


6(2n+l) 





-6(2n+l) 





(2-1) 


-4(3rt 2 -2rt-2) 


4(3n 2 -2n-2) 


-4(3rt 2 -2rt-2) 





4(3n 2 -2n-2) 





(2-2) a 


-2(3rt 2 -2rt-l) 


2(3n 2 -2n+l) 





2(3n 2 -2n) 





-2(3n 2 -2n) 


(2-2) b 





4(6n 2 -4n-l) 








4(3n 2 -2n-2) 


-4(3n 2 -2n) 


(2-2) c 


2(3n 2 -2n-2) 


-2(3ri 2 -2n-2) 


2(3n 2 -2n-2) 





-2(3n 2 -2rt-2) 





(2-3) a 





2(4n-l) 








2(2n+l) 


-2(2n-l) 


(2-3) b 


-(2n+l) 


(2n+l) 


-(2n+l) 





(2n+l) 





(2-3) c 


-2n 


(2n-2) 





(2n-l) 





-(2n-l) 


(3-1) 


-6(2n-2) 


6(2n-2) 


-6(2n-2) 





6(2n-2) 





(3-2) a 


(2n-l) 


-(2n+l) 





-2n 





2n 


(3-2) b 


(2n-2) 


-(2n-2) 


(2n-2) 





-(2n-2) 





(3-2) c 





-2(4n-l) 








-2(2n-2) 


4n 


(3-3) a 














-1 


+ 1 


(3-3) b 


1 








1 

2 





1 

2 


(3-3) c 


-5 





5 
2 





5 
2 





Total sub 


8(3n 2 -2n-l) 





5/'(n) 




-3/'(n) 


-3/'(n) 



Table 15: iA c f x i[S'] n >i — I. All contributions are multiplied by 
lS^£^>i(f) n . Here /'(n)=(6„ 2 -4„-f). 
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(I-J)sub 


7°7 fe Ax fc 
Ax 6 


Ar; 

Ax 6 


Ha~f°Ar]~f k Ax k 


Ai s 


tfaA?? 2 
Ax 6 


Ha'Arj 2 
Ax 6 


(1-1) 


16(n 3 -n 2 -2n) 


-16(n 3 -n 2 -2n) 


16(n 3 -n 2 -2n) 





-16(n 3 -n 2 -2rt) 





(1-2) 


-4(n 3 -n 2 -2n) 


4(n 3 -n 2 -2n) 


-4(n 3 -n 2 -2n) 





4(n 3 -n 2 -2n) 





(1-3) 


6(n 2 +n) 


-6(n 2 +n) 


6(n 2 +n) 





-6(n 2 +n) 





(2-1) 


-4(ra 3 -n 2 -2n) 


4(n 3 -n 2 -2n) 


-4(ra 3 -n 2 -2n) 





4(n 3 -n 2 -2ra) 





(2-2) a 


— 2(rt 3 — n 2 — n) 


2(n 3 — n 2 +n) 





2(n 3 -n 2 ) 





-2(n 3 -n 2 ) 


(2-2) b 





A(2n i -2n 2 -n) 








4(n 3 -n 2 -2n) 


— 4(n 3 — n 2 ) 


(2-2) c 


2(n 3 -n 2 -2n) 


-2(n 3 -rt 2 -2n) 


2(n 3 -n 2 -2n) 





-2(n 3 -n 2 -2n) 





(2-3) a 





2(2n 2 -n) 








2(n 2 +n) 


-2(n 2 -n) 


(2-3) b 


— (n 2 +n) 


(n 2 +n) 


-(n 2 +n) 





(n 2 +n) 





(2-3) c 


-n 2 


(n 2 -2n) 





(n 2 -n) 





-(n 2 -rt) 


(3-1) 


-6(n 2 -2n) 


6(n 2 -2n) 


-6(n 2 -2n) 





6(n 2 -2n) 





(3-2) a 


(n 2 ~n) 


-(n 2 +n) 





-n 2 





n 2 


(3-2) b 


(n 2 ~2n) 


-(n 2 ~2n) 


(n 2 -2n) 





-(n 2 -2n) 





(3-2) c 





-2{2n 2 -n) 








-2(n 2 -2n) 


2ra 2 


(3-3) a 














— n 


+n 


(3-3) b 


n 








-in 





in 


(3-3) c 


— 5n 





-§« 





h 





Total sub 


8(n 3 -n 2 -n) 





5/(n) 


/(«) 


-3/(n) 


-3/(n) 



Table 16: iA c f x i[S'] n >i — I'. All contributions are multiplied by 

^£^>l(f) n in(f). Here /(n)=2n 3 -2n 2 -fn. 
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The same pattern here happened again! We should keep it in mind that this 
pattern might maintain after summing up each individual contribution. We 
also separate temporal terms with spatial ones for our conventional choice 
of counternterms and summarized the results in Tables 15, 16, 17 and 18 
before many summations are performed. Table 16 (Table 18) is the partner 
of Table 15 (Table 17) with the extra logarithm, In |. From the coefficient 
of each individual term at the bottom of each table one might already notice 
that the pattern we mentioned above still exists. The benefit for postponing 
the infinite summation for each individual contraction not only because it 
is a less complicated procedure but also because the pattern serves us one 
consistent check whether or not we have done the computation correctly for 
such a long and complicated calculation. 

We can easily read off the contribution from each distinctive term. For 
example, the total coefficient of 7 \^ Xk from Table 15 and Table 16 is, 



K*mH g J 8(3n 2_ 2n _ 1) + 8(n 3_ n 2_ n) ln( | ; 



64tt 4 



{ y_ y 



(140) 



After collecting all of terms from the four tables we obtain, 



K 2 mH ^ , y ( 7°7 fc Ax k 



64tH 



E 



-4' 



Ax ( 



8(3n 2 -2n-l) + 8(n 3 -n 2 -n) ln(|) 



+ 



,Ha-y°Arry k Ax k Ha'-y Ar]j k Ax k 



Ax 6 



Ax 6 



HaArj 2 Ha'A-q 2 



+ (2n 3 -2n 2 -^n) ln(|) 
H 2 a 2 ^ k Ax k 



Ax 6 
Ha Ha' 



Ax 6 



(6n 2 — 4n — ) 



Ax 4 Ax 4 



(6n 2 +2n) + (2n 3 +n 2 )ln(^) 



Ax 4 
H 2 aa'^ k Ax k 



Ax 4 



-(9n 2 + 12n+ Jj) + (-3n 3 - 6n 2 - \n+ \) ln(|; 

;-3n 2 -2n+- , ) + (-n 3 -n 2 +-W- : ) ln(^) 
v 4 4 4 4 



(141) 



One might notice that the infinite series could be summed easily using the 
following identities, 



n=l - 1 



Y 



Y 



Y 



n=l 



J2n 2 Y n 



n=l 



Y{Y + 1) 

(i-y)3 



(i-Yy 

Y(Y 2 + 4Y + 1) 



; £n 3 F« 



n=l 



[l-YY 



(142) 
(143) 
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(I-J)sub 


Ha 

Ax 4 


Ha' 
Ax 4 


H 2 a 2 -y°~t k Ax k 
Ax 4 


H 2 aa'~/°y k Ax k 


H 2 a 2 A V 
Ax 4 


H 2 a 2 Ha> 
Ax 2 


(1-1) 


-8(3n 2 -l) 





-4(3n 2 +4n-l) 





4(3n 2 +4n-l) 


2(3n 2 +6n+2) 


(1-2) 


2(3n 2 -l) 





(3n 2 +4n-l) 





-(3rt 2 +4n-l) 


-±(3n 2 +6n+2) 


(1-3) 


-3(2n+l) 





-§(2n+3) 





§(2n+3) 


§(2n+3) 


(2-1) 


2(3n 2 -l) 





(3n 2 +4n-l) 





-(3n 2 +4n-l) 


-i(3n 2 +6n+2) 


(2-2) a 


-(3n 2 +2n) 


-2(3n 2 +n) 


-i(3rt 2 +4n-l) 


-(3n 2 +2n) 


-|(3n 2 +4n-l) 


-±(3n 2 +6n+2) 


(2-2) b 


-2(2n+l) 


-2(6n 2 +2n) 








-2(3n 2 +4n-l) 


-(3n 2 +6n+2) 


(2-2) c 


-(3n 2 -l) 





-i(3n 2 +4n-l) 





i(3n 2 +4n-l) 


i(3n 2 +6n+2) 


(2-3) a 





-(4n+l) 








-(2n+3) 


-±(2n+3) 


(2-3) b 


5(2n+l) 





±(2n+3) 





-i(2n+3) 


-±(2n+3) 


(2-3) c 


-i(2n+l) 


-(2n+i) 


-±(2n+3) 


-±(2n+l) 


-±(2n+3) 


-±(2n+3) 


(3-1) 


6(n-|) 


3 
2 


3n 


3 
2 


-3n 


-|(2n+l) 


(3-2) a 


l(4n+l) 


(2n+|) 




i(4n+l) 




|(2n+l) 


(3-2) b 


-(n-l) 


l 

4 




l 

4 


h 


|(2n+l) 


(3-2) c 


2 


(4n+l) 








2n 


|(2n+l) 


(3-3) a 


-1 


1 








l 

2 


l 

4 


(3-3) b 





l 
2 


l 

8 


1 

4 


1 
8 


1 

16 


(3-3) c 


5 
2 





5 
8 





5 
8 


5 

16 


Total sub 


-6(3rt 2 +n) 


-6(3n 2 +n) 


-(9n 2 +12n+|) 


(-3n 2 -2n+|) 









Table 17: iA c{ x ^S^m — II. All contributions are multiplied by 

K 2 mH v^oo ( M.\ n 
64vr 4 2^n>l\i) ■ 
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(I-J)sub 


Ha 

Ax 4 


Ha' 
Ax 4 


Ax 4 


H- 2 aa' 7 °7 fc Ai fc 
Ax 4 


H 2 a 2 A v 
Ax 4 


H 2 a 2 Ha' 
Ax 2 


(1-1) 


-8(n 3 -n) 





-4(n 3 +2n 2 -n-2) 





4(n 3 + 2n 2 -n-2) 


2(ra 3 +3n 2 + 2n) 


(1-2) 


2(n 3 -n) 





(n 3 +2rt 2 -n-2) 





-(n 3 +2n 2 -n-2) 


-j(n 3 +3n 2 +2n) 


(1-3) 


-3(n 2 +n) 





-§(n 2 +3n+2) 





|(n 2 +3n+2) 


|(n 2 +3n+2) 


(2-1) 


2(n 3 -n) 





(n 3 +2n 2 -n-2) 





-(n 3 +2n 2 -n-2) 


-j(n 3 +3n 2 +2n) 


(2-2) a 


— (n 3 +n 2 ) 


-(2n 3 +n 2 ) 


- 1 (n 3 +2n 2 -n-2) 


-(ra 3 +n 2 ) 


- 1 (n 3 +2n 2 -n-2) 


-±(n 3 +3n 2 +2n) 


(2-2) b 


-2(n 2 +n) 


-2(2n 3 +n 2 ) 








-2(n 3 +2n 2 -n-2) 


-(n 3 +3n 2 + 2n) 


(2-2) c 


-(n 3 -n) 





- 3 j(n 3 +2n 2 -n-2) 





|(n 3 +2n 2 -n-2) 


|(n 3 +3n 2 +2n) 


(2-3) a 





-(2n 2 +n) 








-(n 2 +3n+2) 


- ;j(n 2 +3n+2) 


(2-3)b 


\{n 2 +n) 





i(n 2 +3n+2) 





-±(n 2 +3n+2) 


- |(n 2 +3n+2) 


(2-3) c 


-±(n 2 +n) 


-(n 2 + ±n) 


-±(n 2 +3n+2) 


-±(n 2 +n) 


-±(n 2 +3n+2) 


-|(n 2 +3n+2) 


(3-1) 


(3n 2 -fn) 


-§„ 


|(n 2 -l) 


fn 


-f(n 2 -l) 


- 3 (n 2 +n) 


(3-2) a 


±(2n 2 +n) 


(n 2 + ±n) 


i(« 2 -l) 


±(2n 2 +n) 


\{n 2 -l) 


±(n 2 +n) 


(3-2) b 


-(±n 2 -fn) 




-K« 2 -i) 




i(« 2 -l) 


±(n 2 +n) 


(3-2) c 


2n 


(2n 2 +n) 








(„ 2 -l) 


±(n 2 +n) 


(3-3) a 


-("+§) 


(n+i) 








\{n+l) 


\{n+l) 


(3-3) b 


i 

8 






|(2n+l) 


\{n+l) 




(3-3) c 




5 
8 


f("+l) 


5 
8 


-f(n+l) 




Total sub 


-3(2n 3 +n 2 ) 


-3(3n 3 +n 2 ) 


(_ 3n 3_ 6n 2_3 n+ 9 ) 


(-n 3 -« 2 +|T^|) 









Table 18: iA cf x i [•S'jn^i — II'. All contributions are multiplied by 

R 2 mH y^oo fjA™ WiM 
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Here Y stands for |. After the summation, the total contribution from 
iA ci i[S] n >i is, 



— i 



-icfnl 



K 2 mH j 7VAx fc 



64tt 4 



Ax 6 



(~Y 2 +7Y) | (Y*+6Y*-Y) 



(1-YY 



+ 



Ha~f Arrf k Ax k Ha''y Arrf k Ax k HaArj 2 Ha'Ar/ 2 
5 : — - 1 1 — 3—: — - — 3- 



Ax G 
3r 3 +22F 2 -3F 



Arc 6 



2(i-r) 4 

H 2 a 2 ^ k Ax k 



in(y) 



+ 



Ha 



Ax 6 
Ha'' 



Ax 6 



{l-YY 

-3r 3 +26r 2 +r 
2(i-r) 3 



Ax 4 Ax 4 



-12(r 2 + 2F) -3(Y 3 +8Y 2 +3Y) ' 



(i-y) J 



(i-y)< 



Ax 4 
H 2 aa'^ k Ax k 
+ Ax 4 



-3(F 2 -6Y+29) -3(3r 4 -12F 3 +23F 2 +10F) , 

+" ~ — ln(y) 



4(1-F) 3 



4(i-r) 4 



(5y 3 -14F 2 -15F) (-3F 4 +14F 3 -35F 2 ) , 

H 77: — 777-; m(r ) 



4(1-F) 3 



4(1-F) 4 



(144) 



5.2 Sub-Leading Contributions from i5A 



A 



In this subsection we work out the contribution from substituting the residual 
A- type part of the graviton propagator in Table 5, 



a/3 pa 



(x; x') 



VapVafi + VaaVpfi' 



D-3 



Va?V 



pa 



i5A A (x;x') . (145) 



As with the conformal contributions of the previous section we first make the 
requisite contractions and then act the derivatives. The result of this first 
step is summarized in Table 19 and Table 20. We have sometimes broken 
the result for a single vertex pair into as many as five terms because the 
three different tensors in (145) can make distinct contributions, and because 
distinct contributions also come from breaking up factors of 7 a J /?M . These 
distinct contributions are labeled by subscripts a, b, c, etc. 

The next step is to act the derivatives and it is of course necessary to 
have an expression for i5Aa(x; x'). From (66) one can infer, 



i5Aa(x] x 

h 2 r 



D 



1) 



(aa') 2 -it 



167T? "f — 2 



Ax D ~ 4 



+ 



H°- 2 T(D-l) 



H 



(47T) 



(4tt)# r(f) 

D - 2 ^l T(n + D-l) ^yy 1 



7T 



-7rcot(— Dj + ln(aa') 



T(n+f • 



2 n =l 



n T(n+ 



n- 



+ 2 T(n + 2) M 



(!)' 



'-+2\ 



. (146) 
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I 


J 


sub 


iV^x) i[S](x; x>) %V P j a {x') [ afj T£] iSA A (x; x') 


1 


4 




— pz|) iK 2 am$i [S] (x; x')SA A (x; x') 


2 


4 




(d-3) iK 2 am^i [S] (x; x')i5A A (x; x') 


3 


4 


a 


— 2{d^3)^ k2 am l° doiSA A (x; x')i[S](x; x') 


3 


4 


b 


— 2?D-l) in 2 am ffi8A A (x ; x'^jS^x; a;') 


4 


1 




g"Jjm 2 am<9;{i[5](a;; x')7^A A (x; x')} 


4 


2 




^_^iK 2 amdk{i [S] (x; x')^kibA A (x; x')} 


4 


3 


a 


— 2(d^3) i k2 am ^ [S] ( x i ^')7°^o^Aa(x; x') 


4 


3 


b 


-gEl^ami^Kx; x')^A A (x; x') 



Table 19: Contractions from the iSA A part of the graviton propagator-I 



In D = 4 the most singular contributions to (99) have the form, iSA A / Ax . 
Because the infinite series terms in (146) go like positive powers of Ax 2 
these terms make integrable contributions to the quantum-corrected Dirac 
equation (1). We can therefore take D = 4 for those terms, at which point 
all the infinite series terms drop. Hence it is only necessary to keep the first 
line of (146) and that is all we shall ever use. 

The generic contraction from Table 19 which only consists of one deriva- 
tive acting on a propagator, the order m contributions must be and could 
only be from the most singular part of the fermion propagator. In reducing 
these terms the following derivatives occur many times, 

^<5A A (x;x') = _^! r r(^ + l)(aa') 2 ^^^ = _^A A (x;x'),(147) 
8-7TT v 2 ' Ax u A 

mx -, x >) = 5r(f +1 )( aa f-«{^-^L} 

H D - 2 V{D-1) 



2 ■ 

o/.c-a / /x H 2 ( D a, /x2 d[ Art a'H 

& i5A A {x]x) = — u-r( — + l)(aa') 2 ~~ ' 



2Dn , r(Dl oH, (148) 



87r f ^2 ^ ' I Ax^- 2 2Ax^- 4 
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I 


J 


sub 


tVf(x) i[S](x; x>) iVj prT (x') [ Q/J 7£] ^A A (x; x') 


1 


1 






1 


2 




(D ^ 3) K 2 a fc ^[5](x;a; / )7fc^A A (a;;x / )} 


1 


3 


a 


2 ( £ -1 2 \ ^ [S] (x; x')7°d' i5A A (x; x') 


1 


3 


b 


^r7 2 l k 2 &i \S] far; x') &i5A A (x: x') 


2 


1 




1 n 2 d' M\S](x- x')^iSAa(x- x')\ 


2 


2 


a 


. r ^ ft L J V ' / t ^rLV " /J 


2 


2 


b 


+ |/t 2 (947fc^[<S'](^; x / )7fci<5A j4 (x; x')} 


2 


2 


c 


2(d-3) K 2 d k $i{S}(x; x')~f k i5A A (x; x')} 


2 


3 


a 


■?tn ^K 2 $i[S](x; x')<y°d' i5A A (x; x') 


2 


3 


b 


_i K 2 q > i\s]( x - x 'h (k de)i5A A (x; x') 


2 


3 


c 


— „ , * Q .K 2 Bi\S](x: x f )Bi5AA(x: x r ) 


3 


1 


a 


|(^i)« 2 9;{ 7 9o^A A (x;x / )i[-5](x;x')7' 1 } 

^ \ j^j — r* \ ' / L J \ ' / ' J 


3 


1 


b 


i^ K ^iSA A (x; x')i[S](x; x')^} 


3 


2 


a 


2(T } ^K 2 d k {^doi5A A (x] x')i[S](x; x')^ k } 


3 


2 


b 


Ain nsK 2 d k $i5A A (x; x')i[S](x; x')^ k } 


3 


2 


c 


\K 2 ${i[S]{x-x')$i5A A {x-x')} 


3 


2 


d 


\n 2 d k {^i[S\{x\ x')-f e d k i5A A (x; x')} 


3 


3 


a 


-li^^iSKx; x'h°d d' t5A A (x; x') 


3 


3 


b 


j^K 2 l°i[S](x; x')d $i5A A (x; x>) 


3 


3 


c 


-§J^)^[S](x- x')d kl °d' i5A A (x- x') 


3 


3 


d 


iI{d-z) k2 ™[S\{x] x')d k $i5A A (x] x') 


3 


3 


e 


-±K 2 -f k i[S](x; x')-f k V 2 i5A A (x; x') 



Table 20: Contractions from the iSA A part of the graviton propagator-II 
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(I - J)sub 


16^ 2 ( D ~ 4 ) 1 ' 


lJ ^ T ma8 A (x — x') 


(1-4) 


32 i/°- 4 
(D-3) 


-121n(a) 


(2-4) 


n-4 (O-l) 
^ (D-3) 2 


3 ln(a) 


(3 - 4) a 


n-4 (O-l) 
^ 2(D-3) 2 


|ln(a) 


(3 - 4) b 


n-4(0-l)(0-2) 
^ 2(D-3) 2 


3 ln(a) 


(4-1) 





9 + 61n(-g) 


(4-2) 








(4 - 3) a 


n-4 (D-l) 
r 2(D-3) 2 


|ln(a) 


(4 - 3) b 


n-4(0-l)(0-2) 
^ 2(D-3) 2 


31n(a) 


total 


32 H°- 4 n-4D(D-l) 
v^F (D-3) ^ (S-3) 2 


9 + 61n(g) 



Table 21: The local terms from i5A A x i[S] cf . 



We also make use of following identities to simplify the contributions, 

Ar] 2 _ 1 2 1 11 

A"^ 2 ~ 4(£)-2)(£)-3) °Ai 2 ^ ~~ 2(D-2) Ax 2D ^ ' (15 °^ 

A? 1^1 (D-l) 1 n , n 

+ n/r. o\ A~3TO » ( 151 ) 



Arr 2 ^ 2 4(D-2)(D-3) Ax 2 ^ 2( J D-2)Ax 2 ^ 

7 °A?n fc Ax fc _ -1 oo a L n ^ 

A* 2 ^ 2 ~4(D-2)(D-3) 7 ^Ax 2 ^' (i52j 

Ar ? 1 a 1 7 fc Aa; fc -1 5 1 



Ax 21 ^ 2{D-3) v Ax 21 ^ ' Ax 2L ^ 2{D-3y Ax 21 ^ 

Note that the contraction (1-1) produces a delta function through (103) and 
picks up one divergent and one finite local term which do not possesses any 
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dimension-dependent powers of Ax. We tabulate this kind of the result in 
Table 21. One might already notice that the first two terms of (146) 12 and 
the final two terms of (148) and (149) tend to cancel in D = 4. We indeed 
encounter this kind of cancelation entirely for the leading divergent terms in 
contractions (4-1) and (4-2) and hence we present the remaining finite results 
in Table 22. For the rest of the contractions, the sum of the leading divergent 
contributions do not vanish after we apply (150), (151), (152), (153), (127), 
(128), (129), (110) and (111). Therefore we give the local terms in Table 21 
and the finite nonlocal terms in Table 22. After collecting all contributions 
from Table 21 and Table 22 we obtain, 



MdAcf 



\X\X = 



ik 2 H 2 T(^ 



ma 



(D-4) 



32 H 



^D(D-l) 



x5 D (x-x')l 



ik 2 H 2 



X 



ln(/x 2 Ax 2 



Ax 2 



16tt 2 



+ 



9 + 61n(g) 



ma5 



x- 



V^(D-3) 

,. k 2 H 2 
-x') + 



32tt 4 



(D-3) 2 
mal [6<9 2 -2V 5 



/9 nl ,H 2 
( - + 2 n( — 

V2 V 4/x : 



:))V 2 



Ax 2 



(154) 



The generic contractions from Table 20, which are comprised of two 
derivatives acting on the propagators, the order m contributions could either 
come from the flat spacetime mass term or from the infinite series expansion 
of the fermion propagator. 

We first deal with the contributions from the flat spacetime mass term of 
the fermion propagator. At D = 4 these terms have a dimension ^jf^t which 
are not integrable when we substitute back to the quantum-corrected Dirac 
equation and working on an arbitrary D is required. The contractions (1-1), 
(1-2), (2-1) and (2-2) tend to cancel through (146). It turns out that the sum 
of divergent terms vanish and only left the finite terms. Hence we present 
the results for D = 4 at the end. Here we showed (2-1) as an example, 



-K 



-d' , 
(D-3) ^ 

Ax 2D - 



h[S]YiSA A }- 



k 2 h 2 r(f )r(f+i) 



32tt d 
K 2 H D-2 r ( jD _D 



(D- 



2 D+l n D p_3) 



-3) 
-2 



mad 1 ,, 



[aa 



A2- 



_D 

2 



k 2 H 2 pP~Sna 
32n D (D-3) 



(D-4) (D-4) 



(D-4) 

1 , H 2 
— In — - 

2 4/i 2 



f ln(aa') 



\(D-4) 



Ax D 



Ax 2 



\n(fi 2 Ax 2 
Ax 2 



2 7TCOt(^) = T 



(D-4) 
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(I-J)sub 


a 2 MfiAx) 2 
° Ax 2 




V72ln(MAx) 2 
Ax 2 


U Ax 2 




VV-2 

Aa; 2 


(1-4) 




















(2-4) 


3 
4 


-1 


-1 





-2- ln (|J) 


-2-ln(|J) 


(3-4). 


3 
8 








3 
4 


3 
4 





(3-4) b 


3 
4 











1 

2 


i 

2 


(4-1) 


3 

















(4-2) 





1 


-1 









(4-3). 


3 
8 








3 
4 


3 
4 





(4-3) b 


3 
4 











1 
2 


i 

2 


total 


6 





-2 


3 
2 


1 

2 


-l-21n(g) 



Table 22: The non-local terms from ISAa x i[5 r ] c f - All contributions are 
multiplied by ^3 ma. 



(I -J) 


a 2 In(pAx) 2 
° Ax 2 




V/2 M^Aa;) 2 
V Ax 2 


Ffn/) ln (^) 2 




(1-1) 


-3 








3 


3 


(1-2) 





-1 


1 





-1 


(2-1) 





1 


1 








(2-2) 








1 

2 








total 


-3 





5 
2 


3 


2 



Table 23: iSAa x i^fm. All contributions are multiplied by ^^ma. 
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K 2 H 2 

32tt 4 



ma< 



3 
2 



-ln- 



Aft 2 



Ax 2 



^ ln(n 2 Ax 2 



(155) 



The most singular terms in the contractions (2— 3) a , (3—1) and (3—2) are 
integrable after exacting out derivatives 13 . Because the contraction (3-3) 
involves various double derivatives directly acting on i5Aa we would give 
their expressions here, 



d Q df Q i8A A -- 



# 2 r(f + l) f -1 



87rf(aa')T- 2 lA^- 



Ax D 2 { ' 



H 2 aa'An 2 1 H 2 aa' 



d k d iSA A - 



<9 fc <9oio"A A : 



#T(f + 1 
87r"? (aa' 



-(L>-2)A77Ax fc 1 
Ax D 2 { 



2Ax D ~\ 

HaAxk 



•4) 



H 2 T> 



D 



Ax D-2 

Ha'Axk I 



d k dii5A A - 



— i 

8ir 2 (aa 

# 2 n 



(D-2)A^ 1 

Ax D 2 { ' Ax D ~ 2 



D 



{D-2)Ax k A Xl 



5 



ki 



87r 2 (aa') 2 



Ax D 



Ax D - 



(156) 
',(157) 
(158) 
(159) 



Note that the terms with a factor of (D — A) in the equations (156)-(158) 
have dimensionality, which is either A i- 4 or A g_ 3 ■ When they combine with 
if-S'jfm whose dimensionality is A „„_ 2 , one can see that those contributions 



Az D -2, 

are integrable in four dimensions and hence could only give the contributions 
of the order (D — A). Therefore we can drop the terms we mentioned above 
from (156) to (158) when we compute the contraction (3-3). In addition, the 
contractions (3— 3)b and (3— 3) c are both finite in D = A dimensions after per- 
forming the partial integration 14 . The finite contributions mentioned above 
are displayed in Table 23 and Table 24 and they can be read off immediately, 



K 



32vr 4 



-3<9 2 + -V 2 + 3Had + 2Harf? 



In /i 2 Ax 2 
Ax 2 



+ 



tt2 

Sln^d 2 



7 



+ (5 + § h 5) v ' + (l +31n S) a * + ( 1+2h ^)^ 



An 



Aii 2 



A^i 2 " A 

-Ll 

Ax 2 ) 



. (160) 



13 Some less singular terms in the contractions (2— 3) a , (3— l) a and (3— 2) a get canceled 
through the property of (148) and (149). 

14 "Partial integration" is not standard usage in Physics. Here we mean extracting the 
derivatives outside the quantum-corrected Dirac equation. 
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n — i) , 


° Ax 2 


7 °o<7 Ax 2 


V Ax 2 


Hnf).. 1 


nu >i v Ax 2 


(1-1) 


2 6 1 %? 








2 +,3m 4 /1 2 


2 +,31n 4 At 2 


(1-2) 





_3_T i£i 

2 n V 2 


2 +1I V2 





2 "V 2 


(2-1) 





2 +1I V 2 


1 + ^ 








(2-2) 








4 2 4/j 2 








(2-3). 





1 

4 











(3-1). 


3 
2 


3 
2 


3 
2 








(3-l)b 





-1 


-1 








(3-2) 





1 

2 


3 
4 








(3 — 3) b +c 





1 1 

4 4 











total 


-31ng 


7 
4 


2 + 2 m 4 M 2 




l + 2\n$ 



Table 24: i5A^ x i^fm- All contributions are multiplied by |^-ma. 
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(I— J) sub 


Ar, 2 


Ax* 


1 


-y°Arfy k Ax k 


HaArj 


Ha-Y°j k Ax k 


Ax 2D-2 


Ax 2D-2 




Ax 2D ~ 2 


Ax 2D ^ 


Ax 2D^t 


(1-3). 


-(D-l) 








(D-l) 


(D-l) 
(D-2) 





(l-3)b 





-(D-2) 





-(D-2) 





-(D-2) 


(D-3) 


(D-3) 


(D-3) 


(2— 3)b+ c 





(D-l)(D-2) 2 
4(D-3) 














(3-3) a 


(D-l)(D-2) 
2(D-3) 





(D-l) 
2(D-3) 











(3— 3)d+ c 





(D-2) 2 (D-5) 
8(D-3) 


-(D-l)(D-2)(D-5) 
8(D-3) 











total 


(D-l)(3D-8) 
2(D-3) 


3(D-2) 2 (D-5) 
8(D-3) 


-(D-l) 2 (D^i) 
8(D-3) 


-1 

(D-3) 


(D-l) 
(D-2) 


-(D-2) 
(D-3) 



Table 25: The contribution from i5A^ x i[S] im which are not integrable in 
four dimensions. Note that all contributions are multiplied by the factor 
t^r(f-l)r(f+l)ma(aa') 2 ^. 



The rest of the contractions which require further renormalization are 
summarized in Table 25. We then apply the same formalism to partially 
integrate, extract the local divergences and take D = 4 for the remaining, 
integrable and ultraviolet finite nonlocal terms. The sub-total from Table 25 
could be obtained, 



8(D-2)(D-3) 2 (D-4) 



iK 2 H 2 3ma 4 K 2 H 2 ma 
+ — n — - — lnao (x— x ) + 



3 2 \n/j, 2 Ax 2 



64tt 2 2 



32tt 4 32 



d 



+ 



Ax 2 



Ax 2 



(161) 



Combining (160) and (161) gives, 



-ddAfm 



, i\ in 2 H 2 D ,/^(D-l)(Z> 3 -9D 2 +20D-4) 



64ttt v 2 



8(D-2)(D-3) 2 (D-A) 



. n , A iK 2 H 2 3ma, r4 . . K 2 H 2 ma 
mad (x — x)-\ Inao^x— x ) + 



64vr 2 2 



327T 4 



_^9 2 +-V 2 +3Hado 
32 2 
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T T 

1 J sub 




Ar/ 3 
Ax 6 


7°A ) ? 2 7 fe Az fc 


Ar; 
Ax* 


■y°f k Ax k 


Ax 6 


Ax 4 


/ -i -i \ 
(1-1) 


1 








24(3n 2 +4n+l) 


-24(3n 2 +4n+l) 


(1-1) 


lnY 








24(n 3 +2n 2 +n) 


-24(n 3 +2n+n) 


(1-2) 


1 











8(3n 2 -l) 


(1-2) 


lnY 











8(n 3 -n) 


(2-1) 


1 


-16(3n 2 -6n+2) 


16(3n 2 -6n+2) 


-8(6rt 2 -6rt+l) 


8(3n 2 -2n) 


(2-1) 


lnY 


-16(n 3 -3n 2 +2n) 


16(n 3 -3n 2 +2n) 


-8(2n 3 -3n 2 +n) 


8(n 3 -n 2 ) 


(2-2) a 


1 


-2(3n 2 -6n+2) 


2(3n 2 -6n+2) 


-(6n 2 -6n+l) 


(6n 2 -2n-l) 


(2-2) a 


lnY 


-2(n 3 -3n 2 +2n) 


2(n 3 -3n 2 +2n) 


-(2n 3 -3n 2 +n) 


(2n 3 -n 2 -n) 


(2-2) b 


1 


-6(3n 2 -6n+2) 


-2(3ri 2 -6ri+2) 


-3(6ri 2 -6n+l) 


-(6n 2 -2n-l) 


(2-2) b 


lnY 


-6(n 3 -3n 2 +2n) 


-2(n 3 -3n 2 +2rt) 


-3(2n 3 -3n 2 +n) 


-(2n 3 -n 2 -n) 


(2-2) c 


1 


4(3n 2 -6n+2) 


-4(3n 2 -6n+2) 


2(6n 2 -6n+l) 


-2(6n 2 -2rt-l) 


(2-2) c 


lnY 


-4(n 3 -3n 2 +2n) 


-4(n 3 -3n 2 +2n) 


2(2n 3 -3n 2 +n) 


-2(2n 3 -ra 2 -rt) 


total 


1 


-20(3n 2 -6n+2) 


12(3n 2 -6rt+2) 


2(6n 2 +78n+7) 


-6(6n 2 +18n+5) 


total 


lnY 


-20(n 3 -3n 2 +2n) 


12(n 3 -3n 2 +2ri) 


2(2n 3 +39n 2 +7n) 


-6(2n 3 +9n 2 +5n) 



Table 26: i5A A x i[)S'] n >o — IA. The factor 
multiplies all contributions. Here Y 
factors for the each individual row. 



2 H 2 mHaa' H 2 Ax 2 



2 6 7T 4 



Eoo -yr. 
n=0 1 



|; InF and 1 are the multiplicative 



+2//a 7 V 



ln(/x 2 Aa; 2 



Ax 2 

(4 



+ 



+ f?+3 In^)/Ja9 + (2 + 2 ln^)ffo 7 V 



2 
1 



Ax 2 



(162) 



The final class is comprised of terms in which comes from the infinite series 
expansion. Theses contributions are integrable in D = 4 and do not require 
any further renormalizations so we could set D = 4 right at the beginning 
without alternating the final result. Therefore we shall apply the infinite 
series expansion of the fermion propagator for D = 4 to this calculation. It is 
actually the same expression as (137). In addition, one can also use (138) and 
all related identities in Appendix A but remember that the series is summed 
up from n = instead of n — 1. The residual A- type graviton propagator 
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I-Jsub 




ffoAij 2 


ffa'A, 2 


Ha-/ A v -/ k Ax k 


Ha'j°Ari^ k Ax k 


Ha 
Ax 2 


Ax* 




Ax* 


Ax* 


(1-1) 


1 


24(3n 2 +4n+l) 





-24(3n 2 +4n+l) 





12(3n 2 +6n+2) 


(1-1) 


lnY 


24(n 3 +2n 2 +n) 





-24(n 3 +2n 2 +n) 





12(ra 3 +3n 2 +2n) 


(1-2) 


1 








8(3n 2 -l) 








(1-2) 


lnY 








8(n 3 -n) 








(2-1) 


1 


4(3n 2 -l) 


-4(3n 2 -l) 


-4(3n 2 -l) 


4(3n 2 -l) 


2(6n 2 +6n+l) 


(2-1) 


lnY 


4(n 3 -n) 


— 4(n 3 — n) 


-4(n 3 -n) 


4(n 3 -n) 


2(2n 3 +3n 2 +n) 


(2-2)a 


1 


(3n 2 -l) 











i(6n 2 +6n+l) 


(2"2) a 


lnY 


(n 3 — n) 











±(2n 3 +3n 2 +n) 


(2-2) b 


1 


3(3n 2 -l) 











§(6n 2 +6rt+l) 


(2-2) b 


lnY 


3(n 3 — n) 











|(2n 3 +3n 2 +n) 


(2-2) c 


1 


-2(3n 2 -l) 











-(6rt 2 +6n+l) 


(2-2) c 


lnY 


-2(n 3 -n) 











-(2n 3 +3n 2 +n) 


total 


1 


6(15n 2 +16n+3) 


-4(3n 2 -l) 


-4(15n 2 +24n+7) 


4(3n 2 -l) 


9(6n 2 +10n+3) 


total 


lnY 


6(5n 3 +8n 2 +3n) 


-4(n 3 -n) 


-4(5n 3 +12n 2 +7n) 


4(n 3 -n) 


9(2n 3 +5n 2 +3n) 



Table 27: i5A A x i[S] n > - IB. The factor mH 2 aa ' (in ^f^j En=o Yn 
multiplies all contributions. Here Y — |; lnY and 1 are the multiplicative 
factors for the each individual row. 
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I-Jsub 




Ha' 
Ax 2 


H 2 a 2 Ari 
Ax 2 


H 2 a 2 1 % k Ax k 
Ax 2 


H 2 aa' 1 ° 1 k Ax k 
Ax 2 


H 3 a 2 a' 


(1-1) 


1 





-6(3rt 2 +6rt+2) 


6(3n 2 +6n+2) 





-3(3n 2 +8rt+5) 


(1-1) 


lnY 





-6(n 3 +3n 2 +2n) 


6(n 3 +3n 2 +2n) 





-3(n 3 +4n 2 +5rM-2) 


(1-2) 


1 








-2(3n 2 46n+2) 








(1-2) 


lnY 








-2(n 3 +3n 2 +2n) 








(2-1) 


1 


-2(6n 2 +6n+l) 








-2(3n 2 +4n+l) 





(2-1) 


lnY 


-2(2n 3 +3n 2 +n) 








-2(n 3 +2n 2 +n) 





(2-2)a 


1 

















(2-2). 


lnY 

















(2-2) b 


1 

















(2-2) b 


lnY 

















(2-2) c 


1 

















(2-2) c 


lnY 

















total 


1 


-2(6n 2 +6n+l) 


-6(3n 2 +6n+2) 


4(3n 2 +6n+2) 


-2(3n 2 +4n+l) 


-3(3n 2 +8n+5) 


total 


lnY 


-2(2n 3 +3n 2 +n) 


-6(n 3 +3n 2 +2n) 


4(n 3 +3n 2 +2n) 


-2(n 3 +2n 2 +n) 


-3(n 3 +4n 2 +5n+2) 



Table 28: i5A A x i[S} n > - IC. The factor ™ 2 g 2 mH 2 aa ' (in ^f^j En=o Yn 
multiplies all contributions. Here Y — |; YaY and 1 are the multiplicative 
factors for the each individual row. 
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I-J 




Ar, 3 
Ax 6 


7 A ) 7 2 7 fe Ax fc 
Ax 6 


Ar/ 
Ax 4 


Ax 4 


(1-1) 


1 








12(3n 2 +8n+3) 


-12(3n 2 +8n+3) 


(1-1) 


lnY 








12(n 3 +4n 2 +3n) 


— 12(n 3 +4n 2 +3n) 


(1-2) 


1 











4(3n 2 +4n+l) 


(1-2) 


lnY 











4(n 3 +2n 2 +ra) 


(1-3) 


1 








-6(2n+l) 


8(2n+l) 


(1-3) 


lnY 








-6(n 2 +n) 


8(n 2 +n) 


(2-1) 


1 


-8(3n 2 -2n) 


8(3n 2 -2n) 


-4(6n 2 +2n) 


4(3n 2 +2n+l) 


(2-1) 


lnY 


-8(n 3 -n 2 ) 


8(n 3 -n 2 ) 


-4(2ra 3 +n 2 ) 


4(n 3 +n 2 +rt) 


(2-2) 


1 


-2(3n 2 -2n) 


-2(3n 2 -2ra) 


-(6n 2 +2n-3) 


-(6n 2 +6n-l) 


(2-2) 


lnY 


-2(n 3 -n 2 ) 


-2(rt 3 -n 2 ) 


-(2n 3 +n 2 -3n) 


-(2n 3 +3n 2 -ra) 


(2-3) 


1 


10(2n-l) 


-6(2n-l) 


2(10rt-2) 


-2(2n+l) 


(2-3) 


lnY 


10(n 2 -n) 


— 6(n 2 — n) 


2(5n 2 ~2ra) 


-2(n 2 +n) 


(3-1) 


1 


-8(2n-2) 


8(2n-2) 


2(4n-l) 


-2(8n-3) 


(3-1) 


lnY 


-8(n 2 -2n) 


8(n 2 -2n) 


2(2n 2 -ra) 


-2(4n 2 -3n) 


(3-2) 


1 


-2(2n-2) 


-2(2n-2) 


-(4n-l) 


(4n-3) 


(3-2) 


lnY 


-2(n 2 -2n) 


-2(n 2 -2n) 


-(2n 2 -n) 


(2n 2 -3n) 


(3-3) 


1 


-5 


3 


7 

2 


l 

2 


(3-3) 


lnY 


— 5n 


3n 


-h 


-in 


total 


1 


-5(6n 2 -4n-l) 


3(6n 2 -4n-l) 


(6n 2 +98rt+f ) 


-(18n 2 +78n+f ) 


total 


lnY 


-5(2n 3 -2n 2 -n) 


3(2n 3 -2n 2 -n) 


(2n 3 +49n 2 +f n) 


-(6n 3 +39n 2 +^n) 



Table 29: x i[S] n > - IIA. The factor mH 2 aa ' £^ =0 ^™ multiplies all 
contributions. Here Y — |; lnY and 1 are the multiplicative factors for the 
each individual row. 
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I-J 




HaAr/ 2 
Ai 4 


Ha'Ar, 2 
Ai 4 


Ha-f°Ari~f k Ax k 
Ax 4 


Ha'~i Q Ari'y k Ax k 
Ax* 


Ha 
Ax 2 


(1-1) 


1 


12(3n 2 +8n+3) 





-12(3n 2 +8n+3) 





6(3n 2 +6n+2) 


(1-1) 


lnY 


12(n 3 +4n 2 +3n) 





-12(n 3 +4n 2 +3n) 





6(n +3n +2n) 


(1-2) 


1 








4(3n 2 +4n+l) 








(1-2) 


lnY 








4(n 3 +2n 2 +n) 








(1-3) 


1 


-6(2ra+l) 





8(2n+l) 








(1-3) 


lnY 


-6(n 2 +n) 





8(ra 2 +n) 








(2-1) 


1 


2(3n 2 +4n+l) 


-2(3n 2 -l) 


-2(3n 2 +4n+l) 


2(3n 2 -l) 


(6n 2 +14n+5) 


(2-1) 


lnY 


2(n 3 +2n 2 +n) 


-2(n 3 -n) 


-2(n 3 +2n 2 +n) 


2(n 3 -n) 


(2n 3 +7n 2 +5n) 


(2-2) 


1 


(3n 2 +4n+l) 











i(6n 2 +14n+5) 


(2-2) 


lnY 


(n 3 +2n 2 +n) 











i(2n 3 +7n 2 +5n) 


(2-3) 


1 


-3(2n+l) 





2(2n+l) 





-3(2n+l) 


(2-3) 


lnY 


-3(n 2 +n) 





2(n 2 +n) 





-3(n 2 +n) 


(3-1) 


1 


4n 


-2(2n+l) 


-2(2n) 


2(2n+l) 


-(2n+3) 


(3-1) 


lnY 


(2n 2 -2) 


-2(rt 2 +n) 


-2(rt 2 -l) 


2(n 2 +n) 


-(n 2 +3n+2) 


(3-2) 


1 


6n 





2(2n) 





6n+f 


(3-2) 


lnY 


3(n 2 -l) 





2(n 2 -l) 





3n 2 + fn+f 


(3-3) 


1 


5 
2 











7 

4 


(3-3) 


lnY 


f(n+l) 













total 


1 


5(9n 2 +20«+^) 


-2(3n 2 +2n) 


-2(15n 2 +34n+12) 


2(3n 2 +2rt) 


27n 2 +55r)+x 


total 


lnY 


5(3n 3 +10n 2 +f n-i) 


-2(n 3 +n 2 ) 


-2(5n 3 +17n 2 +12n) 


2(n 3 +n 2 ) 


9n 3 +f n 2 +fr^§ 



Table 30: i5A A x i[S] n > - IIB. The factor mH 2 aa ' £~ =0 yn multiplies all 
contributions. Here Y — |; lnY and 1 are the multiplicative factors for the 
each individual row. 
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I-J 




Ha 1 
Ax 2 


H 2 a 2 Ari 
Ax 2 


H 2 a 2 j% k Ax k 
Ax 2 


H 2 aa'~f°~f k Ax k 
Ax 2 


H 3 a 2 a' 


(1-1) 


1 





-3(3n 2 +10n+8) 


3(3n 2 +10n+8) 





-§(3n 2 48r!rf5) 


(1-1) 


lnY 





-3(n 3 +5n 2 -|-8r!+4) 


3(n 3 +5n 2 -|-877rf4) 





-|(n 3 +4n 2 45n+2) 


(1-2) 


1 








-(3n 2 +10«+8) 








(1-2) 


lnY 








-(n 3 +5n 2 +87j+4) 








(1-3) 


1 





§(2n+3) 


-2(2n+3) 








(1-3) 


lnY 





§(n 2 +3«f2) 


-2(rt 2 +3n+2) 








(2-1) 


1 


-(6n 2 +6«+l) 








-(3n 2 +4n+l) 





(2-1) 


lnY 


-(2n 3 +3n 2 +n) 








-(n 3 +2n 2 +n) 





(2-2) 


1 

















(2-2) 


lnY 

















(2-3) 


1 

















(2-3) 


lnY 

















(3-1) 


1 


(2ra-2) 








2(2n+2) 





(3-1) 


lnY 


(n 2 -2n-3) 








2(ra 2 +2rH-l) 





(3-2) 


1 

















(3-2) 


lnY 

















(3-3) 


1 

















(3-3) 


lnY 

















total 


1 


-(6rt 2 +4«+3) 


-3(3n 2 +9rM-^ ) 


2(3n 2 +8n+5) 


-(3n 2 -3) 


-§(3n 2 4«r!+5) 


total 


lnY 


-(2n 3 +2rt 2 +3nf3) 


_ 3(n 3 + | n 2 + M r)+ 3 ) 


2(n 3 +4n 2 45ref2) 


-(n 3 -3n-2) 


-|(n 3 +4n 2 +5n+2) 



Table 31: i5A A x i[S] n > - IIC. The factor %f, mH f J2n=o yn multiplies all 
contributions. Here Y — |; lnY and 1 are the multiplicative factors for the 
each individual row. 
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and its various derivatives in four dimensions occurred very frequently, 



iSA A (x; x) 



-H 2 r H 2 Ax\ h 

H—^ + 2 



8tt 2 



dii5A A (x; x') = —dii5AA(x;x') ■ 

d i5A A (x;x') = —d' i5A A (x;x') 
—H 2 1 

d d' i5A A (x;x') 



Air 2 



+ 



-H 2 Axi 

An 2 Ax 2 
H 2 At] 

~ An 2 Ax 2 
2A V 2 - 



Ax 2 Ax 4 



d dii5AA(x;x') = —d' diiSAA(x;x') 



H 2 -2Ar]Axi 
4^ Ax 1 



didii5A A (x;x') 



H 2 r-S, 2A Xi Ax, 



An 2 



Ax 2 



+ 



Ax 4 



(163) 
(164) 
(165) 
(166) 
(167) 
(168) 



We also make use of the following identities to facilitate our computation 
more effectively, 



/ k 

11 ll 



i k , Yi°i» 



2 7 °, (169) 
7 y 7 VAx fc = i Q i k Ax k , 7 fc 7 i 7°7 / Aa; fc = -?>i°i k Ax k , (170) 
7V7V'Aa7 = 3A V + i°i k Ax k , (171) 
7°7 M 7°7 iy A^Ax !/ = Ax 2 - 2 7 °A^Aa; /1 , (172) 
l k ri°l u Ax k Ax l ,Ax^-i i k Ax k Ax 2 -2Ar n k Ax k YAx lx , (173) 
l°l u l°l k Ar]Ax u Ax k = ArjAx 2 - ^Arf^Axk , (174) 

(175) 

— ) and 



1 k 1 u 1 ° 1 l Ax k Ax l Ax v = ArjAx -rf^ k Ax k Ax . 
Note that any derivatives acting on i5A A would eliminate ln(- 



that the exceptions are the generic contractions (1-1), (1-2), (2-1) and (2-2). 
We list theses results in Table 26, 27 and 28. The rest of the contributions 
without In h2 ^ x2 are summarized in the Table 29, 30 and 31. From these 
tables one can see that the derivative of the coefficient with In | is the coeffi- 
cient without In |. Based on the characteristic of (137) we should not be too 
surprised at this peculiar pattern occurring here again as in the previous sub- 
section. The final result for — i[£ dAn0 ](a;; x') could be computed using (142) 
and (143) and then adding the lowest order constant to each distinct con- 
tribution because the summation here starts from n = rather than n = 1. 
Finally we tabulate the lengthy results in Table 32 and Table 33. 
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ln H 2 Ax 2 
4 


1 


Ar/ 3 
Ax 6 


-20 A (Y) - 20 fo(Y) lnY 


_5o (y) _ 5g 2 (y) i n y 


7°Ar/ 2 7 fc Ax fc 
Ax 6 




<JyH-< J + og 2 \i ) in -i 


Ar/ 
Ax* 


9 f„('V"| -|- 9 f, (Y) In Y 


nofy'i _i_ n . (Y\ ln y 
i/3 V 1 ) \ i/4 \ 1 1 in 1 


7°7 fc Ax fc 
Ax 4 




nr-(V\ nJV\ ln V 
95 { 1 ) <?6 V 1 ) m 1 


HaA-n 2 
Ax 4 


vj7\i ) -T VJ8[I ) 111 I 


y7\ 1 1 ' °y8\ 1 ) HI -1 


Ha 1 A-q 2 
Ax 4 


—4fn(Y) —4f\n(Y) lnY 


— 2onfyi — 2omfyi iny 


Haj°Ar)~f k Ax k 
Ax 4 


A f (V\ A f (V\ In V 




Ha'-y°Arn k Ax k 
Ax 4 


4/ 9 (y) + 4/ 10 (y)iny 


2g 9 (Y) + 2g l0 (Y) lnY 


Ha 

Ax 2 


9/i 3 (y) + 9/ 14 (y) lnY 


9is(Y)+gu(Y) lnY 


Ha' 

Ax 2 


-2/ 15 (Y) - 2/ 16 (Y) lnY 


-gi 5 (Y) - g 16 (Y)\nY 


H 2 a 2 A v 
Ax 2 


-6/ 17 (y) - 6/ 18 (y) lnY 


-3g 17 (Y)-3g 18 (Y) lnY 


H 2 a 2 j°-y k Ax k 
Ax 2 


4/ 17 (y) + 4/ 18 (y)i n y 


2g 19 (Y) + 2g 20 (Y) lnY 


H'^a'^^Axk 
Ax 2 


-2/ 19 (Y) - 2f 20 (Y) lnY 


-92i(Y) - g 22 (Y) lnY 


H 3 a 2 a' 


-3/21 (Y) - 3/ 22 (y) lnY 


-|^ 9 (y) - f<? 20 (y) iny 



Table 32: The total result for i5A A xi[S] n > . The factor m ^, Qa ' multiplies 
all contributions. Here Y — |; ln h2 ^- x2 and 1 are the multiplicative factors 
for the each individual column. The various functions fi{Y) and <7j(Y) are 
presented in Table 33 
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£ f~\.r\ 

fi(Y) 




9i(Y) 




hi Y ) 


Y(2Y 2 +5Y-1) . 




-Y(Y 2 -12Y-1) i 


(l-y) 3 1 z 


(i-y) 3 


f2{Y) 


6Y 3 

(1-Y) 4 


92{Y) 


Y(3Y 2 +10Y-1) 

(i-y)« 


hv) 


Y(7Y 2 -86Y+91) , 7 

; 1 L -I— / 


93{Y) 


y(49y 2 -282y+257) , 49 


(i-y) 3 1 ' 


2(i-y) 3 1 2 


n{Y) 


-6y(5y 2 +y-8) 


94{Y) 


y(-45Y 2 -82Y+151) 


(i-y)4 


2(l-y) 4 




Y(5Y 2 -22Y+29) . r 




y(37y 2 -194y+229) 37 


(i-y) 3 1 ° 


2(1-Y) 3 1 2 


£ r\ 


2Y(-Y 2 -Y+8) 
(1-Y)* 


<?e(30 


y(-29y 2 -26y+i27) 
2(i-y) 4 


fi(Y) 


Y(3Y 2 -7Y+34) , Q 




y('13y 2 -48Y+71~) l^i 


(i-y) 3 1 6 


2(1-Y) 3 1 2 




2Y(7Y+8) 
(l-y)4 


0s(5O 


y(y 3 -4y 2 +y+38) i 
2(i-y) 4 2 


h(Y) 


Y(-Y 2 + 5Y+2) -, 

(i-y) 3 1 




y(y+5) 

(i-y) 3 


M Y ) 


6y 2 
(i-y)4 


9w(X) 


2y(2y+i) 
(i-y) 4 


fn(Y) 


y(7y 2 -23y+46) 7 


9n(Y) 


y(i2y 2 -43y+6i) 10 


(i-y) 3 1 ' 


(i-y) 3 1 iz 


M Y ) 


6Y(Y+4) 
(I-Y)* 


9i2{Y) 


-2y(2y-i7) 
(i-y) 4 


M Y ) 


y(3y 2 -ioy+i9) . Q 


9n{Y) 


y(79y 2 -27oy+407) . 79 


(1-Y)3 1 <5 


4(i-y)3 1 4 


M Y ) 


2y(y+5) 


9u{Y) 


Y(-5Y 3 +20Y 2 -29Y+230) , 5 


(i-y)4 


4(l_y)4 1 4 


M Y ) 


y(y 2 -2y+i3) . -, 


9i5[Y) 


y(3y 2 -4y+i3) , Q 


(i-y) 3 1 1 


(l-y)3 1 ° 


M Y ) 


6y(y+i) 


9w{Y) 


y(-3y 3 +i2y 2 -7y+io) . Q 


(i-y)4 


(l_y)4 1 <J 


M Y ) 


y(2y 2 -7y+n) , 

— — — 1— / 


9i7{y) 


y(i3y 2 -38y+37) . 13 


(l_y)3 1 z 


2(i-y) 3 1 2 


MY) 


6y 
(i-y) 4 


9is(Y) 


3Y(-Y 3 +4Y 2 -6Y+5) , Q 

(i-y) 4 1 


MY) 


y(y 2 -3y+8) , •, 


9w(Y) 


y(5y 2 -i5y+i6) , c 


(i-y) 3 1 1 


(l-y)3 1 


My) 


2y(y+2) 


92o(Y) 


-2Y(Y 3 -4Y 2 +6Y-6) , 


(i-y) 4 


(i-y) 4 1 Z 


MY) 


y(5y 2 -i5y+i6) , r 


92i(Y) 


-3Y 2 (y-3) Q 


(1-Y)3 1 d 


(i-y) 3 6 


MY) 


-2y(y 3 -4y 2 +6y-6) , 


922{Y) 


2Y(Y 3 -4Y 2 +8Y-2) r> 


(l-y)4 1 z 


(i-y) 4 A 



Table 33: The coefficient functions for the table 32 
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5.3 Sub-Leading Contributions from i8As 

In this subsection we work out the contribution from substituting the residual 
B- type part of the graviton propagator in Table 5, 



a/3 A pa 



As in the two previous sub-sections we first make the requisite contractions 
and then act the derivatives. The result of this first step is summarized in 
Table 34. We have sometimes broken the result for a single vertex pair into 
parts because the four different tensors in (176) can make distinct contribu- 
tions, and because distinct contributions also come from breaking up factors 
of 7° ? J /3m . These distinct contributions are labeled by subscripts a, b, c, etc. 

25Ab(x; x') is the residual of the B-type propagator (67) after the confor- 
mal contribution has been subtracted, 



lfor* Ai D - J ( 4jr )f r(f) 

JjD-2 oo fr/„x5 



- 2 ^ f T(n+f ) /yxn-f +2 T(n+D-2) ,y,n] 
lk\nn+2)yV " r(n+f) ^ j ^ 



As was the case for the iSAa^x') contributions considered in the previous 
sub-section, this diagram is not sufficiently singular for the infinite series 
terms from i8A B (x; x') to make a nonzero contribution in the D = 4 limit. 
Unlike idA^ix; x'), even the n = terms of iSA B (x; x') vanish for D = 4. Most 
of the contractions involve at least one derivative of iSA B , 

d t i5A B (x; x') = - H2V ^f (D-4)(aa') 2 -% = -%5A B (x; x% (178) 

167T 2 

d 0i SA B (x-,x') = ^^(D-4)(aar-H^-^^\^ (™) 



d Q & i5A B (x;x') 




# i5A B (x;x>) = ^^(D-4)(aaT^{-^-7^-A i (180) 

167T 2 



ff 2 r(| ) (p-4) 

167rf (aa')^ 2 l Aa;D A ^ 2 

+£>[CD-4) 2 ] , (181) 



d d k iSA B (x;x') = -d' d k 8A B (x;x') 



64 



I 


J 


sub 


iVj^Cx) i[S](x; x>) iVF(x') [ Q/3 T*] ^A B (x; *0 


2 


1 







2 


2 


a 


-|/t 2 ^{ 7 (°9 fc )i[ 1 S](x; x')7 fe i<5A B (x; x')} 


2 


2 


b 


-|K 2 9 fc {7 (0 9 fc) i[51(x; x')7°^A B (x; a/)} 


2 


3 


a 


-|K 2 7fc9Di[S'](x; x')7 fc d' Q i5A B (x; x') 


2 


3 


b 


\K 2 ^°d' i5A B (x] x') ^[^(x; x') 7 fc 


2 


3 


c 


-\K 2 ^ k d k i5A B (x\ x') d i[S](x] x') 7 ° 


2 


3 


d 


§K 2 7°d fe i[S](x; x')7° d k iSA B {x; x') 


3 


1 







3 


2 


a 


lK 2 d' {>y k i[S}(x; x')^ k d i5A B (x; x')} 


3 


2 


b 


lK 2 -f k d k {i[S]{x; x')7° d iSA B {x; x')} 


3 


2 


c 


-|« 2 7 °9;{z[5](x; x') 7 fc ^A B (x; x')} 


3 


2 


d 


-lK 2 d k {~?°i{S}(x; x')7° d fe ;<5A B (x; x')} 


3 


3 


a 


-±n 2 'y k i[S\(x; x')<y k d d' i5A B (x; x') 


3 


3 


b 


^ 2 7°i[S](x; xV d k & i6k B {x; x') 


3 


3 


c 


-^^[S] (x; x')7° dodkiS^Bix; x') 


3 


3 


d 


±K 2 7°i[S](x; x')7 V 2 2(5A B (x; x') 



Table 34: Contractions from the iSA B part of the graviton propagator. The 
contributions from (l-l)-(l-4), (4-l)-(4-4), (2-4) and (3-4) are zero. 
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pre — factor 


K 2 H 2 rfD\TtD i ^ ma 


K 2 ffM r(D~2)r(fi) 




2-£>+4 7r -D r(— ) 


(2-2) ai 




-(z>-i)[^+ J H"a9 ]^^2 


(2-2) a2 






(2-2) bl 






(2-2) b2 


-(^2)^2 1 
2(D-3) V Ax 2 "- 6 


V7 2 1 



Table 35: The contributions for the contraction (2-2) from iSA B x ifS'jfm. 



ff 2 r(|) (g-4)(g-2) i -ArjAxk 
16ttt (aa')^" 2 I Ax D 



■+0[(D-4) 2 ], (182) 



o PA , M H 2 T(D)2 (D-A) \(D-2)Ax 2 (D-l)} . . 
V^A B (^) = ^ ( ^ |L_±_-L_i|. ( 18 3) 

The fact that the first line of i5AB(x;x') and its various derivatives are 
all of the order (D — A) means that they would only contribute when they 
are multiplied by a divergence. Note that the contractions consisting of mass 
interaction vertices all vanish through (176), so no order m contributions 
come from i5A^ x ifiSjcf. The potential non-zero, order m contributions 
could either come from the flat spacetime mass term or from the infinite 
series expansion of the fermion propagator. Remember that the only term in 
^[S^fm behaves like Ax d-2 and that the most singular one in i[>S] n >o goes like 
Ax d-3 - In addition, the generic contractions in Table 34 are comprised of two 
derivatives. Therefore one can count that the dimensionality of most singular 
contribution from ifS*]^ is . 20-4 whereas the one from i[S'] n >o is » 20-5 • 

The 

former is logarithmically divergent in D = A before performing the partial 
integration, so one still needs to keep arbitrary D for the computation; the 
latter is entirely integrable in D = A so that one could compute this part in 
four dimensions and the result turns out to be zero owning to the cancelation 
of the first two series of i5Ab(x;x') and owning to (D — A) factor from its 
various derivatives. Therefore the only class we need to work out in this 
sub-suction is iSA B (x;x') x i[S] im (x; x'). 

We take special care of the contraction (2-2) because it is the only contrac- 
tion in Table 34 which derivatives might have a chance not to act upon iSAs- 
We also break up ^°d k ^ into ^■y°d k and |7 fc <9° for each sub-contraction and 
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I Jsub 


contributions 


(2-3) a 




(2-3) b 




(2-3) c 


2(D-3) ' a 0? Ax 2D-6 


(2-3) d 




(3-2) a 


(D-4)(D-1) pp. 1 
(£>-3) U Ax 2D ~ 6 


(3-2) b 


(r) _8)7 Cofe^ 


(3-2) c 


(M„,0a 5 1 

(D-3) > ^9 Ax -2D-6 


(3-2) d 


(£^4)^72 1 

(D-3) V Ax 2D - 6 


(3-3) a 


-(0-1)^2 1 
4(D-3) a Ax 2D - 6 


(3-3) b 


-(M.,Oa a 1 
4(D-3) ' °W A:E 2i;- 6 


(3-3) c 


ig|)7%fe^ 


(3-3) d 


r-(^-l)fl2 , (H V721 1 
L 4(£>-3) W 1 4(£>-3) V J Ax 2D - 6 


total 


\2(D-3)° 1 ( D 4 )L(D-3) a oJ 4(D-3) V J A^"« 



Table 36: i5A B x i[S , ] fm terms. All contributions are multiplied by 
^mar(f)r(f-l)(aa') 2 -T. 
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tabulate the results in Table 35. Theses expressions are integrable in four di- 
mensions and the contributions from the left column cancel out exactly with 
that from the right column in D = A. Table 36 gives the rest of our results 
for the most singular contributions, those in which all derivatives act upon 
the conformal coordinate separation. There is no net contribution when one 
or more of the derivatives acts upon a scale factor. Those expressions are 
also integrable in D = A dimensions, at which point we can take D = 4 and 
the result vanishes on account of the overall factor of (D — A) or (D — 4) 2 . 
Finally we read off the net contribution from Table 36 and take D — A, 



-udBfm 



. k 2 H 2 3 n2 1 k 2 H 2 3 Q2 1 , 10/| , 

( X ; x') = -r—rma< -d 2 -— } = ——-ma< -d 2 -— } . (184) 
v ' ; 2 10 vr 4 |2 Ax 2 | 2% 4 I 8 Ax 2 | v ; 



5.4 Sub-Leading Contributions from i5A 



c 



The point of this subsection is to compute the contribution from replacing 
the graviton propagator in Table 5 by its residual C-type part, 



apA pa 



^2 



Vapfjp* SZtyfjpr + fjgpSffi / D 3 \ q q q q 

+ — wz-o — + ta J 5 *Wo 



i5A 



c- 



(D-2)(D-3) ~ D-2 ^\D-2r aPpa 

(185) 

As in the previous sub-sections we first make the requisite contractions and 
then act the derivatives. The result of this first step is summarized in Table 37 
and Table 38. We have sometimes broken the result for a single vertex 
pair into parts because the four different tensors in (185) can make distinct 
contributions, and because distinct contributions also come from breaking 
up factors of 7 a J /?M . These distinct contributions are labeled by subscripts 
a, b, c, etc. 

Here i5A c (x;x') is the residual of the C-type propagator (68) after the 
conformal contribution has been subtracted, 

... , „ H 2 ( D , ( D ^{aa') 2 -% H D - 2 T{D-3) 

i6A c (x;x') = w ( 3)r( 1) \ ' - + w ,ns 

1 ^ 167TT ^2 ) V 2 J Ax (47T)# T(f) 



H 



(4vr)f n tll V 2 I T(n+2) ^ ^ ' I>+f ) W J 



(186) 



As with the contributions from i5As{x\ x') considered in the previous sub- 
section, the only way i5Ac(x; x') can give a nonzero contribution in D = A 
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I 


J 


sub 


tV^(x) i[S](x; x>) iVr(x>) UTfJ i5A c ( X] a;') 


1 


4 




(g-3) 2 (D-2) K2 ^ am ffi [S] ( x 'i x')i5A c (x; x') 


2 


4 


a 


(d-2) K ' 2 'i a ' m l doi[S](x; x')i5Ac(x; x') 


2 


4 


b 


— ( D _J S ( D _ 2 ) f^iam @i[S] (x; x')i5A c (x; x') 


3 


4 


a 


^^^^ia mi d i5A c (x-x')i[S}(x-x') 


3 


4 


b 


2(d-3)(d-2) K2 i am $i8Ac(x; x')i[S](x; x') 


4 


1 




(D-3KD-2) ^am^jilSKi; a^^cte a/)} 


4 


2 


a 


— p32j/c 2 iam^{i[<S'](a;; x')7°i<5A c (a;; x')} 


4 


2 


b 


(d-3) 1 (d-2) K 2 iamd k {i[S](x; x')~f k i5A c (x; x')} 


4 


3 


a 


2(d-3)(d-2) fi2wm i I s ] ( x 'i x 'h d iSA c (x; x 1 ) 


4 


3 


b 


2(d-3)(d-2) ft2 ^ Qm ^ [^l x> ) 0SA c (x;x') 



Table 37: Contractions from the i5A c part of the graviton propagator-I. 



dimensions is for it to multiply a singular term. That means only the n = 
term can possibly contribute. Even for the n = term, both derivatives must 
act upon the coordinate separation to make a nonzero contribution in D = 4 
dimensions. 

All the vertex pairs involve one or more derivatives of i5Ac, 
d li 6A c = H2T i 2n 1 \ ^-3)(D-4)(aar^ ^ = -8[r5A c , (187) 



lGvrf V 2 \aa')T-i\ Ax D Ax D ^j 

+0[(D-A)% (190) 

d k d iSA c = -d k d' i5A c 
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J 


sub 


iVf(x) i\S](x: x') iV? a (x') \ a8 T c } iSA c (x: x') 


1 


1 




-7dw-2t« 2 ^,{^[5](^; x')7^a c (x;x')} 


1 


2 


a 


-mb)^ W^](x; ^)7°^A c (x; x')} 


1 


2 


b 


rn sin 2 ^ 2 <9jfc{^[<S](x;x07 fe ^A c (x;x')} 


1 


3 


a 


2<dVd 2)* 2 x ; h°d' t5A c (x; x') 


1 


3 


b 


~27n wn 2^ 2 ^[^K^ x ') ^<5A c (x; x') 


2 


1 


a 


m 1 9^ 2 <9' {7°<9o^](x; x')7^A c (x; x)} 

I — Z ) fit K 1 ^ L J V ' / / w \ 3 /J 


2 


1 


b 


m ^^d'ABilS}^: x'^iSAcix: x)\ 

(D — 6)(D — 2) // Lr L J \ ' ) i ^ ' /J 


2 


2 


a 


2 %- 3 2 \ « 2 ^{7°d i[S](x; x')7°^A c (x; x')} 

z( D — z ) U l ' u L J \ ' /I \ ? /J 


2 


2 


b 


9m 1 ^ /t 2 <9fc{7°<9o^l(x; x')7 fc i<5A c (x; x')} 


2 


2 


c 


i^^S-]^; x')7°^A c (x; x)} 


2 


2 


d 


2m wn 2 r^ 2 ^fc{^['5'](^; x')7^<Ac<(x; x)} 


2 


3 


a 


I(n~9> 2 7°do*[S](x; xV^A^x') 


2 


3 


b 


~4m oT^T^t' 5 ']^; ^0 @iSA c (x; x') 

l± t 1 / ^ 1 ' " L J \ ' / ' \ 1 / 


2 


3 


c 


~Z77^Wn 9^ 2 ^[^(x; x'^c^Ac^x; x') 


2 


3 


d 


aid Wn 9^ 2 ^[S] ( x ; x ^A c (x; x') 


3 


1 


a 


- 2775S)7lh) ^{7% *<5A c (x; (x; x') 7 M } 


3 


1 


b 


2(D-3)(D-2) « 2 ^{^A c (x; x / )«[.S](x; x')7 At } 


3 


2 


a 


-i(S5|) k25 o{7 ^A c (x;x')^5](x;x')7 } 


3 


2 


b 


-4(lb) K25 o{^A c (x; x')i[5](x; x') 7 °} 


3 


2 


c 


-4(g- g 3Kg-2)« 2 ^W ^A c (x; x')i[S] (x; x') 7 fc } 


3 


2 


d 


4(d-3)(d-2) ^ 2 ^{^Ac(x; xO^-SKx; x')7 A; } 


3 


3 


a 


3^^)«: 2 7^^](x; x')7°9 ^5A c (x; x') 


3 


3 


b 


8(d-3)(d-2) k2 7°^[^] (x; x')-f k d d k t5A c (x; x') 


3 


3 


c 


8(D- g 3Kg-2)"V«[g](^; x')7°4^A c (x; x') 


3 


3 


d 


8(d-3)(d-2) K2 7 fe? [<S] (»; x / )7'<9fe<9i^A c (x; x') 



Table 38: Contractions from the iSAc part of the graviton propagator-II. 
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~ 16^ ( ^" 3) (aaO-- 2 A^ + ° [{D - A) ^ (191) 

dkdlt6Ac = ^^^- 3) Ja^\^~ D (192) 

Note that i5Ac and its various derivatives have the same behaviors as %5Ab- 
The propagator itself tends to cancel in D = 4 dimensions and its various 
derivatives all carry (D — 4) factor. This means that they could give the 
non-zero contributions only when they are multiplied by the singular terms. 
For the generic contraction I in Table 37, the only order m contribution 
must come from the conformal part of the fermion propagator whereas the 
generic contraction II in Table 38 the order m contribution could be either 
from the flat spacetime mass term or from the infinite series expansion of 
the fermion propagator. The fortunate thing is that the contribution from 
the infinite series expansion of the fermion propagator vanishes. The most 
singular terms from this particular contribution have dimensionality Ax 2ds , 
which are integrable in D = 4. Therefore they completely vanish in D = 4 
dimensions owing to the behaviors of the residual part of the C-type graviton 
propagator we mentioned above. 

We first work out the contributions from Table 37, which contain a loga- 
rithmic divergence. The contraction (1—4) and (2— 4) a are simple owning to 
the special property of the conformal part of the fermion propagator (103) 
and hence they only pick up the constant part of iSA c , 

t 1 = 2 %# r ( #) (D-^D-sr ^'-^' (193) 
P-4- = S^(sbf^*-*>- (194) 

The same procedure is employed to make the expressions integrable in D = 4. 
We summarized the contractions (4—1), (4— 2) a and (4— 2) b in Table 39 and 
listed the rest of this category in Table 40. The two series in Table 39 from 
each contraction cancel out with each other precisely in four dimensions. 
Finally the total summation from (193), (194), Table 39 and Table 40 is 
quite simple in D = 4, 

*0 = ^f™^ - + ' (195) 
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I Jsub 



2 5 7r D i l 2 ^ i ^2 L ) (D-3) 



ma 



k H 



2 i>fl 7r D 



T(D-3)ma 



(4-1) 



^ 2 



(D-2)(Z>-3) L/ Ax 2D -« 



5 2 



(D-2) 2 (_D-3) L/ Ax D - 2 



(4-2), 



i 



2(D-2) 



'OA^-2 



(4-2) t 



-i 



i 



2(D-2)(D-3) Ax 2D ~6 



(D-2) 2 (D-3) Ax D ~ 2 



1 



Table 39: Contributions from Table 37 for i5A c x i[S] cf (x; x')-l. 



I Jsub 


Q2 1 


(2-4) b 


-1 
4 


(3-4) a 


-1 

8 


(3-4) b 


-1 

8 


(4-3) a 


-1 

8 


(4-3) b 


-1 

8 



Table 40: Contributions from Table 37 for i5Ac x i[S] c f(x; x')-ll. All the 

(3-4?)(£>-l) 



terms are multiplied by fg^r(f )r(f -1) (D _ 2 2 )2(D _3)2 



ma 
• 



72 



I Jsub 



k 2 H 2 t2(D 

2^TT D 



r 2 (f-i)(3-f 7 



1-11 



-(D-3) 



rd 2 



2Ha~f° 



1 



(D-2) (D-3) '"J Ax 2D "« 



2[-a 2 +gg 7 V] 1_ 

(D-2) (D-3) Ax D 



(1-2). 



-i 



2(D-3) 



Ha 
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So 



1 



A:E 2D-6 



(D-2) 
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i 



L2(D-3) 2 



+ 



1 



(D-2) (D-3) Ax 2D - 



-[9+Harf] g l_ 

(D-2) (D-3) P Ax D - 



(2-1). 



b(D-3) ' a ° (D-2) J 



1 



(D-2) J r Ax 2 



(D-2) "Ai D - 2 



(2-l)b 



i 



2(D-3) 2 



-1 



(D-2) (D-3) ! 



Ax 1 



(2-2), 



[-H 2 -^)^ 



Ax°- 



(2-2), 



ti(D^7% + S^ 



2(D-2) J^ Ax 2 ^ 



-[ 7 °i9o+Ha7 ] a 1 

2(D-2) ^Ax"- 2 



(2-2) t 



i 



2(D-2) 



i 



Ax D-2 



(2-2) d 



-l 



i 



4(D-3) 2 v Ai 2D - 



2(D-2)(D-3) v Ax 1 



Table 41: Contributions from Table 38 for i5A c x ^^(a; x')-l. 
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(l~3)a 

(l~3)b 
(2-3), 
(2-3) d 
(3-3) a 



<9 2 



i 

Ax 2 °- (i 

-1 
2(D-3) 

-1 
2(D-3) 



1 

4(D-3) 

if — ) 

8VI3-3/ 



(3-3) d 



i 



8(D-3) 



total 



(D-g) (D-4) 
8(D-3) 4(D-3) 



Table 42: Contribution from Table 38 for i5A c x i[S]f m (x; x')-ll. All the terms 



are multiplied by 



K 2 _ff 2 _r2 / 



r 2 (f-i) 



(3-§)(D-l) 



(^2)(D-3) (aa')«- 
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The final class we need to complete is the contributions from Table 38. 
Very similarly to what happened with i5A#, the contractions from (1-1), (1- 
2), (2-1), (2-2) tend to cancel and we summarized them in Table 41. We 
also tabulate the rest of the contributions which do not vanish in D = 4 
in Table 42. As already explained, terms for which one or more derivative 
acts upon a scale factor make no contribution in D = 4 dimensions, so the 
final nonzero contribution come from the derivatives only acting upon the 
coordinate separation, Ax 2 . The net contributions from Table 41 do vanish 
completely in D = 4 dimensions and the only non-zero contributions of this 
class come from Table 42, 

-i[^ ]M= ^ ma ^j_y (196) 



6 Renormalization 



Except for the finite results from (144), and the results from Table 33-32, 
each of which possesses distinctive expressions, the rest of regulated result 
we have worked so hard to compute derives from summing expressions (98), 
(114), (117), (124), (126), (136), (154), (162), (184), (195) and (196) 

-2 



iK 2 \(3 l -d 2 + P 2 mHd + /3 3 mH^+/3 4 mH 2 a \ 5 D (x -x') + I ( 3 In a - - ) -8 
[a J 167T 2 [ v 8 7 a 

+ (^lna-^)miJ9o+(-^lna + i)m/J7 |+(y kia+^)H 2 ma\s%x-x') 



16 

k 2 r 



16 

3 m 



16 

7 a 27 



/« - I •>///.> fiu — i \ /9a 9 \ ,.,.(].-: > /-J ■) 

64?{ 57* Ms" 32) mff Ml6" 32) m ^ + " M32 



-215 



+ 



32 



+ 



d 2 + H 2 ma 



V 2 + 6Had Q +4Ha-f°0 



\n(fi 2 Ax 2 
Ax 2 



k 2 H< f 9 ff)2 

< — ma o n 

64tt 4 I 16 



-53 
IF 



-ma 



+ ^-ma')-f°d o 0+ (-^ + ln ^)maV 2 + (^ + 6 In ^—)Hma 2 d 
16 7 v 16 4£r 7 v 2 4/i^ 7 



V R R I 4^2/ V 



8 8 



In — J Hmaa 7 



Ax 



2 ' 



(197) 



The various independent constants in (197) are, 
,d-4 r(R_u r_ 2 ( j D_i)] 



>\ — 



r(f-i) 



167rf(£>-3)(£)-4)l (L>-2) 



(198) 
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r(f-i) [ -4(d-i; 

167tt(D-3)(D-4)1 (D-2) 



+ ( J D-2)(6 2 +6 3 



(199) 



D-4 



r(f-i) 



16ttt (£>-3)(£>-4) 



i^_i) + (£,_2)(6 aB +6 3a ) + ^Sd 1 \ (200) 



D-4 



r(f-i) r (p-i) 

16tt# (£>-3)(D-4)\ 2 



+ ( J D-2)(6 4 -6 2 



CD-2)(d 2 + d 3 + d 4 ) 



£>(£>-!) r 



8(D-3) 

(2D — 3) Z) D 



D(L>-2)-l-~(D-4) 



+ 



D-4 



-T{D+1) 



—r^-urp-fj + Qg^j. 



(201) 



Here 6 2 , &2a, &3, &3a, &4, c?ij cfe, d 3 and d 4 are defined in (118) and (125). In 
obtaining these expressions we have always chosen to convert finite, D = 4 
terms with d 2 acting on 1 /Ax 2 , into delta functions, 

' iATi 2 5\x-x') . (202) 



Ax 5 



All such terms have then been included in those which are proportional to 
SHx—x'). 



X 

Fig. 3: Contribution from counterterms. 

The local divergences in this expression are canceled by the BPHZ coun- 
terterms enumerated at the end of section 3. The generic diagram topology 
is depicted in Fig. 3, and the analytic form is, 

4 

_ i r S ctm 

1=1 

m , 
a 

In comparing (197) and (204) it would seem that the simplest choice for the 
coefficients a,- is, 



') = Y J iC Il3 5 D (x-x') , (203) 
i=i 

iK 2 \a 1 —d 2 +a 2 mHdo + a 3 mH^+a 4 H 2 ma}6 D (x-x') . (204) 

l. n > 
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lb 
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Here F 2 
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P-3) 



(205) 
(206) 



Except for the distinctive contributions from (144), Table 33 and Table 32, 
the rest of our final result for the renormalized self-energy is, 
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7 Discussion 

We have used dimensional regularization to compute quantum gravitational 
corrections to the fermion self-energy at one loop order in a locally de Sitter 
background. Our regulated result is (197). Although Dirac + Einstein is 
not perturbatively renormalizable [5] we obtained a finite result (207) by 
absorbing the divergences with BPHZ counterterms. 

For this 1PI function, and at one loop order, only four counterterms are 
necessary. None of them represents redefinitions of terms in the Lagrangian 
of Dirac + Einstein. Two of the required counterterm operators (85) come 
from generally coordinate invariant fermion bilinears of dimension six (73). 
The other two counterterm operators (86) are from other fermion bilinears 
of dimension six (84) which respect the symmetries (44-49) of our de Sitter 
noninvariant gauge (43). 

Although our renormalized result could be changed by altering the finite 
parts of the three BPHZ counterterms, this does not affect its leading behav- 
ior in the far infrared. It is simple to be quantitative about this. Were we to 
make finite shifts A«j in our counterterms (205) the induced change in the 
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renormalized self-energy would be, 



—i 



A£ r 



(x; x') = in 2 l Aa 1 —d 2 +Aa 2 mHd + Aa 3 mH^+a 4 H 2 ma \S 4 (x-x'). 

(208) 



No physical principle seems to fix the Acti so any result that derives from 
their values is arbitrary. This is why BPHZ renormalization does not yield 
a complete theory. However, at late times (which accesses the far infrared 
because all momenta are redshifted by a(t) = e Ht ) the local part of the 
renormalized self-energy (207) is dominated by the large logarithms, 

m 2 { , m 971na — 91na n - 95 met 9 ] rA . ,. . 

3 lna—d 2 + ——mHd + ——mH^+ -^—H^a \ 5%x -x'). (209) 



16tt 2 \ a 16 u 16 ' r 8 

The coefficients of these logarithms are finite and completely fixed by our 
calculation. As long as the shifts A«j are finite, their impact (208) must 
eventually be dwarfed by the large logarithms (209). 

None of this should seem surprising, although it does with disturbing 
regularity. The comparison we have just made is a standard feature of low 
energy effective field theory and has a very old and distinguished pedigree 
[10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 51, 52]. Loops of massless 
particles make finite, nonanalytic contributions which cannot be changed by 
local counterterms and which dominate the far infrared. Further, these effects 
must occur as well, with precisely the same numerical values, in whatever 
fundamental theory ultimately resolves the ultraviolet problem of quantum 
gravity. That is why Feinberg and Sucher got exactly the same long range 
force from the exchange of massless neutrinos using Fermi theory [12, 13] as 
one would get from the Standard Model [13]. 

So we can use (207) reliably in the far infrared. The point of this exercise 
has been to study the effect of breaking conformal invariance with small 
fermion mass. Obtaining (207) completes the first part in that program. 
What remains is to use our result to solve the quantum-corrected Dirac 
equation (1). We shall undertake that in a subsequent paper. However, it 
seems clear that the dominant effect must come from the terms which possess 
large logarithms in local terms and in (144), Table 32 and Table 33 15 . 



15 Without an explicit calculation wc cannot determine whether or not the hi(|) in the 
non-local terms will produce infrared enhancements because the same term occurring in 
[4] fails to give them. 
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As adumbrated in the Introduction, these terms are only enhanced by 
a factor of ln(a) relative to the classical part of the Dirac equation (1). 
That seems roughly the same as the ln(a) enhancement which soft virtual 
gravitons induce on massless fermions [1, 2]. Note that any such effect will 
be independent of assumptions about the existence and couplings of light 
scalars during inflation. 

We have already commented on the importance of the logarithm terms. 
During inflation these infrared logarithms are ubiquitous in quantum cor- 
rections from massless, minimally coupled scalars and gravitons. It is not 
even possible to exclude the possibility that infrared logarithms can contam- 
inate the power spectrum of cosmo logical density perturbations [86, 87, 88]! 
The proportional correction they make in that case must be small because 
the logarithms would only start to grow at horizon crossing, and must cease 
growing when the mode reenters the horizon after inflation. So the largest 
enhancement for a currently observable mode would be ln(a) < 100. This 
must be set against the tiny loop counting parameter of GH 2 < 10~ 10 . 

The more significant corrections would be to modes which are still enor- 
mously super- horizon. These are also down by the constant GH 2 , but the 
time-dependent enhancement factor ln(a) could be arbitrarily big. That is 
what we shall study in our follow-up work. Of course loops of such modes 
can also engender large corrections to effective interactions of low dimen- 
sion. These corrections can grow so large that perturbation theory eventually 
breaks down. Starobinskii has advocated gaining quantitative control over 
this regime by summing the leading infrared logarithms at each order [89]. 
With Yokoyama he has given a complete solution for the case of a minimally 
coupled scalar with arbitrary potential which is a spectator to de Sitter in- 
flation [90]. This powerful non-perturbative technique has been successfully 
generalized to Yukawa theory [33], which showed that the system decays in 
a Big Rip singularity, and to SQED [82], which confirmed the conjecture 
by Davis, Dimopoulos, Prokopec and Tornkvist that super-horizon photons 
acquire mass during inflation. 

The asymptotic late time effect is small in the simple scalar models and in 
SQED for which the leading logarithm expansion have been summed. How- 
ever, the same kind of effect derived from Yukawa is huge. Therefore it is by 
no means clear what might be the outcome for more complicated theories 16 

16 The theory possesses derivative interaction vertices which cannot be avoided by im- 
posing a special gauge as was done in SQED. 
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that also show infrared logarithms such as quantum gravity [70, 71, 72]. An- 
other application of our result (207) is to serve as "data" in checking the 
validity of the new, more general rule [92, 84] for reproducing the leading 
logarithms of massive Dirac + Einstein. This might serve as an important 
intermediate point in the difficult task of generalizing Starobinskn's tech- 
niques to full blown quantum gravity. 

It is well to close with a comment on accuracy. Although parts of this 
computation are quite intricate we have good confidence that (207) is correct 
for two reasons. First, there is the flat space limit of taking H to zero while 
taking the conformal time to be 77 = —e~ Ht /H with t held fixed. This checks 
the leading conformal contributions. Our second reason for confidence is the 
fact that all divergences can be absorbed using just the four counterterms 
we have inferred in section 3 on the basis of symmetry. This was by no 
means the case for individual terms; many separate pieces must be added to 
eliminate other divergences. 
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A The reduced fermion propagator and its 
related identities 

Here we listed some identities we have used for various gamma functions con- 
tracted with the first derivative of the n = part of the fermion propagator 
(122). 
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H 2 a 2 ~f° mH D ~ 3 ,. D . , L> 2 2 n l 

( 7J Zi)A^iJ- W (««')- I r(- + i)r(2--)^ v] 



.(211) 



To facilitate the calculation from the infinite series expansion of the 
fermion propagator for D = 4, we might employ the following identities, 
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The formulae keep the same form for the summation starting from n = 1 as 
long as we are working in D = 4. 
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